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Preface 

Poor  high  angle-of-attack  stability  and  control  characteristics  of 
most  high-performance  aircraft  costs  the  Air  Force,  millions  of  dollars 
annually  through  loss  of  control  accidents.  Even  the  newest  ox  the 
operational  fighter  aircraft  with  their  collection  of  strakes,  slats, 
fences,  boundary  limiters,  sophisticated  flight  control  systems  end 
enormous  vertical  tails,  continue  to  experience  loss  of  control.  The 
report  of  the  Stall/Post-Stall/Spin  Symposium  makes  it  apparent  that 
scientists,  engineers,  and  officials  of  the  Air  Force  and  other  related 
agencies  were  aware  of  the  need  for  an  improved  control  technology  on 
which  to  base  design,  development,  and  testing  of  future  aircraft  and 
to  solve  the  problems  of  current  operational  aircraft.  One  possible 
improvement  is  the  use  of  thrust  generated  control  moments  and  forces 
to  augment  the  aerodynamic  forces  in  the  high  angle-of-attack  regime. 

It  was  my  intention  in  this  study  to  determine  analytically  the  feasi¬ 
bility  of  thrust  control  augmentation.  If  substantially  improved 
handling  and  departure  characteristics  in  the  model  coul  ,  be  demon¬ 
strated,  they  would  serve  as  a  basis  or  incentive  for  f  cure  studies 
of  actual  hardware.  I  was  also  interested  in  improving  the  digital 
computer  program  that  solves  the  non-linear  equations  of  motion  to 
make  it  easier  to  understand  and  use,  simpler  to  modify,  and  to  have 
better  printed  and  plotted  outputs. 

I  want  to  express  my  gratitude  to  my  thesis  advisor,  Lt  Col 
Frederick  F.  Tolle,  Associate  Professor  and  Deputy  Head  of  the  Depart¬ 
ment  of  Aero-Mechanical  Engineering,  for  his  interest,  assistance,  and 
encouragement.  I  wish  also  to  thank  Capt  Donald  C.  F-ckholdt  and  the 
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personnel  of  the  Aircraft  Dynamics  Group,  Air  Force  Flight  Dynamics 
f  '  Laboratory.  I  am  forever  indebted  to  Capt  Eckholdt  for  his  willing 

and  enthusiastic  guidance,  counseling  and  ideas,  and  his  enduring 
patience.  For  my  wife,  there  are  no  words  to  adequately  express  my 
«  appreciation  for  her  patience,  encouragement,  understanding,  and  good 

coffee  during  our  past  four  years  in  AFiT  programs. 


Deas  H.  War ley  III 
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Abstract 

The  feasibility  of  engine  thrust  vectoring  for  lateral  control  of 
aircraft  in  the  high  angla-of-attack  regime  was  investigated  for  an 
airplane  with  F-lll  characteristics.  The  technique  was  found  to  be 
effective  in  increasing  the  angle-of-attack(  at  which  departure  occurs. 
The  method  used  an  effective  dynamic  directional  stability  parameter  to 
account  for  thrust  effect  alteration  of  the  static  lateral  stability 
parameters  Cng  and  Cfcg.  Although  the  effective  Cng  could  not  be 
used  to  predict  departure  in  the  model  studied,  it  was  useful  in  eval¬ 
uating  the  effectiveness  of  the  thrust  vectoring  concepts. 
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ANALYSIS  OF  CONTROL  SURFACE  AUGMENTATION 
IN  HIGH  PERFORMANCE  AIRCRAFT  BY  THRUST  VECTORING 

I.  Introduction 

Most  modern  high-performance,  aircraft  when  flying  at  high  angles- 

of-attack  exhibit  poor  lateral-directional  stability  and  control 

r 

characteristics  which  can  lead  to  inadvertent  departure  or  loss  of 
control.  The  problem  is  to  maintain  control  of  the  aircraft  after  it 
has  lost  aerodynamic  lateral-directional  stability  until  the  classical 
stall  angle-of-attack  has  been  exceeded. 

Background 

Loss  of  control  (departure)  normally  occurs  when  a  pilot  uninten¬ 
tionally  exceeds  the  normal  flight  envelope  boundaries  while  performing 
a  maneuver.  The  situation  is  particularly  probable  and  critical  in 
maneuvers  such  as  landing,  air-to-air  combat,  and  precision  weapons 
delivery  because  of  extreme  pilot  workload  and  stress  level.  Severe 
loss  of  control  may  lead  to  operational  restrictions  which  reduce  the 
mission  effectiveness  of  the  weapons  system. 

Air  Force  safety  records  indicate  tha."  accidents  due  to  loss  of 
control  account  for  almost  25%  of  all  aircraft  losses  and  almost  60% 
of  all  aircrew  fatalities.  The  cost  of  these  losses  is  estimated  to  be 
more  than  40  million  dollars  annually  (Ref  23).  As  of  September  1972, 
at  least  10  of  the  22  F-lll  accidents  (aircraft  destroyed)  were 
attributed  to  loss  of  control.  At  an  estimated  unit  cost  of  20  million 
dollars,  the  price  tag  for  the  F-lll  loss  of  control  accidents  is  a 
staggering  200  million  dollars. 
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In.  the  late  1950’s  (Ref  29)  the  consensus  among  contractors  and 
the  military  was  to  trv  to  solve  the  problem  by  more  f ally  exploring 
the  spin  characteristics  of  each  aircraft,  and  by  teaching  pilots 
proper  spin  recovery  techniques.  Subsequently,  it  has  been  found  that 
most  high-performance  aircraft  have  a  myriad  of  complex  spin  modes, 
several  of  which  are  non-recoverable.  The  loss  of  control  accident 
trend  of  recent  years  and  the  complexity  of  the  spin  problem 'led  to 
the  expected  shift  apparent  at  the  1971  Stall/Post-Stall/Spin  Symposium 
(Ref  5).  Emphasis  is  now  placed  on  either  preventing  the  spin  by  elim¬ 
inating  the  possibility  of  departure  or  controlling  the  spin  by  develop¬ 
ing  and  installing  spin  recovery  devices. 

Despite  advances  in  analysis,  synthesis,  development,  and  testing, 
aircraft  frequently  reach  the  fully  operational  status  before  departure 
problems  are  discovered.  At  that  point,  the  expense  of  a  complete 
research  and  development  program  or  a  large  scale  modification  program 
usually  leads  to  a  less  than  optimal  solution.  The  corrective  action 
has  been  to  train  the  pilot  or  limit  him  by  regulation  in  order  to 
avoid  the  problem-  or  make  minor  configuration  modifications  in  order 
to  mitigate  the  problem. 

Several  examples  of  modifications  to  the  airframe,  control  system, 
or  both  are  described  in  Appendix  E.  Few  of  the  modifications  imple¬ 
mented  to  date  have  appreciably  improved  the  directional  stability 
characteristics  of  the  aircraft.  Additionally,  all  of  the  devices 
discussed  have  at  least  one  of  the  following  drawbacks: 

1.  The  effectiveness  cf  the  device  is  limited  to  a  particular 
flight  regime. 

2.  Margin-of-safety  zones  created  at  the  performance  limits  deny 
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full  use  of  the  maneuver  envelope. 

3.  Pilots  object  strongly  to  any  devices  which  limit  cr  override 
their  control  of  the  aircraft. 

Spin  control  by  using  spin  recovery  devices  is  an  excellent  con¬ 
cept  but  is  limited  by  the  fact  that  there  is  insufficient  reliable 
spin  data  to  evaluate  thoroughly  such  a  system.  Data  acquisition  is 
costly  in  both  time  and  money,  is  complicated  by  the  existence  of 
several  unique  spin  modes  for  each  aircraft  configuration,  and  the 
extensive  flight  testing  required  could  lead  to  further  aircraft  losses. 

These  facts  lead  to  the  conclusion  that  aircraft  must  be  designed 
to  be  directionally  stable  at  angles-of-attack  up  to  and  just  beyond 
the  classical  stall.  Directional  stability  and  control  can  be  achieved 
by  careful  aerodynamic  design  of  the  external  geometry  of  the  airframe 
end  control  surfaces,  by  the  addition  of  nor.conventional  control  sur¬ 
faces  or  forces,  or  by  using  advanced  active  flight  control  systems. 
Directional  stability  through  aerodynamic  design  alone  has  beer,  achieved 
on  only  one  modern  high-performance  aircraft,  the  Northrop  F-5/T-38. 
Since  the  majority  of  high-performance  aircraft  do  not  have  aerodynamic 
directional  stability,  the  only  alternative  that  will  provide  stability 
without  sacrificing  maneuverability  is  to  provide  lateral-directional 
control  by  non- aerodynamic  devices. 

Approach  and  Scope 

The  objective  of  this  study  was  to  investigate  the  feasibility  of 
using  thrust  generated  yawing  moments  and  side  forces  to  augment  the 
aerodynamic  lateral  control  at  high  angles-of-attack.  The  study  was 
limited  to  a  single  aircraft,  the  F-lll,  which  has  directional  stabil¬ 
ity  characteristics  typical  of  many  modern  high-performance  aircraft. 
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A  computer  program  was  developed  (Appendix  D)  to  model  the  air¬ 
craft  using  a  six  degree-of-freedora,  non-linear  formulation  of  the 
equations  of  motion.  It  was  based  on  a  similar  program  developed  at 
the  Air  Force  Flight  Dynamics  Laboratory  (Ref  6) .  The  non-linear 
aerodynamic  data  obtained  from  NASA  was  based  on  wind  tunnel  and  model 
flight  test  data  modified  to  match  the  performance  characteristics  of 
the  F-lll.  While  the  computer  program  car.  be  used  for  spin  analysis, 
the  data  is  valid  only  up  to  departure;  therefore,  this  study  is  limited 
to  investigating  the  application  of  thrust  vectoring  to  prevent  depart¬ 
ure. 

The  approach  is  to  determine  lateral  stability  criteria  as  a 
function  of  angle-of-attack  and  sideslip,  and  to  evaluate  the  influence 
of  thrust  cctrol  augmentation  on  this  criteria.  This  permits  an 
analytic  demonstration  of  the  feasibility  of  thrust  control.  A  dis¬ 
cussion  of  methods  of  implementing  thrust  control  and  a  short  synopsis 
on  the  state-of-the-art  is  presented;  however,  actual  hardware  design 
is  not  included. 
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II.  Development  of  the  Model 


Aerodynamic  Loss  of  Control 

The  boundary  of  the  aircraft  maneuver  envelope  can  be  defined  by 
the  angle-of-attack  that  corresponds  to  either  maximum  lift  or  zero 
dynamic  directional  stability.  Configuration  A  shown  in  Figure  1 
maintains  positive  dynamic  directional  stability  at  all  angles-of- 
attack.  The  aircraft  will  stall  in  the  classical  sense  when  the  angle- 
of-attack  of  maximum  lift  is  exceeded,  and  will  not  experience  uninten¬ 
tional  departure  or  loss  of  control.  In  contrast,  configuration  B  will 
depart  from  controlled  flight  when  the  dynamic  directional  stability 
parameter  (Cn^  ^  )  approaches  zero  prior  to  the  conventional  stall 

point.  Since  this  departure  occurs  at  an  angle-of -attack  that  is  often 
considerably  less  than  that  of  maximum  lift,  the  result  is  the  loss  of 
a  large  portion  of  the  high  angle-of-attack  maneuvering  capability.  It 
is  in  this  angle-of-attack  range  that  the  forces  and  moments  generated 
by  thrust  can  be  used  to  augment  those  of  the  control  surfaces  whose 
effectiveness  has  been  diminished  due  to  adverse  air  flow  character¬ 
istics. 


Model  Equations  of  Motion 

The  system  is  modeled  by  the  boly  axis,  non-linear,  six-degree-of- 
freedom  equations  shown  in  Appendix  B.  The  forces  (X,  Y,  and  Z)  and 
the  moments  (L,  M,  and  H)  are  expanded  as  functions  of  the  aerodynamic 
and  thrust  effects.  Two  general  categories  of  thrust  control  augmenta¬ 
tion  are  considered.  The  first  deals  with  cases  where  the  augmenting 
thrust  system  produces  or  exerts  both  control  moments  and  control 
forces  on  the  aircraft  such  as  that  produced  by  deflection  of  the 
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engine  thrust.  The  second  category  applies  to  cases  where  the  augment¬ 
ing  thrust  system  produces  only  control  moments. 

Engine  Exhaus t  Deflection.  The  model  assumes  that  engine  exhaust 
can  be  deflected  by  scheduling  the  exhaust  nozzles  for  asymmetric 
closing  or  by  flow  separation  induced  by  injection  of  bypass  or  bleed 
air  at  the  nozzles.  Although  the  actual  hardware  is  not  discussed, 
the  technology  exists  and  is  currently  being  applied  to  VTOL  aircraft 
(Ref  3,  8,  and  26).  Engine  thrust  is  assumed  to  be  an  external  force 
acting  on  the  aircraft  at  the  nozzle  position.  The  lengths  x^,  y^,  and 
Z.J,  define  the  location  of  each  nozzle  in  the  body  axis  reference  frame. 
Tp  and  T2  are  respectively  the  left  and  right  engine  thrust  forces. 
Defining  a^,  to  be  the  thrust  vector  angle  with  respect  to  the  x-z 
plane  and  0T  to  be  the  thrust  vector  angle  with  respect  to  the  x-y 
plane,  the  thrust  can  be  expressed  at  any  general  angular  position. 
Maximum  thrust  vector  angles  are  restricted  to  within  the  thrust 
deflection  limits  (15°  maximum)  achievable  by  ai‘-  injection  techniques. 


Resolving  the  forces  and  noments  from  Figure  2, 

X^,  *=■  (^  +  T2)  cos  aT  cos  P>T  '  (1) 
Yt  -  (T1  +  T2)  sin  Sj  (2) 
ZT  -  -  (Tx  +  I,)  sin  ctp  (3) 
Lt  -  (Tx  -  T2)  yT  sin  -  (T1  +  T2)  z T  sin  (4) 
Hj  ■  -  (Tx  +  T2)  xt  sin  aT  +  +  T2)  z T  cos  aT  (5) 
NT  -  -  ai  +  T2)  xt  sin  0T  +  (T:  -  T2)  yT  cos  $T  (6) 


Auxiliary  Thrusters .  Using  a  reaction  control  system  (RCS)  util¬ 
izing  high  pressure  compressor  bleed  air  or  small  rockets  in  coupled 
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FIGURE  2.  Definition  of  Thrust  Auguentation  Position  arid  Angles. 
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pairs  (for  example:  four  nozzles  located  in  opposing  pairs  at  each 
wingtip  and  aligned  parallel  to  the  x  body  axis) ,  control  moments  can 
be  generated  without  the  presence  of  translation  producing  forces. 
Assuming  the  auxiliary  thrusters  act  only  in  the  x-y  body-axis  plane, 
the  result  from  Figure  2  is  the  yawing  moment  equation, 

Nwr  "  (wtx  "  \*  ^  ^ 

Programmed  Equations  of  Motion.  Equations  (1)  through  (7)  are 

combined  with  the  non-linear  equations  developed  in  Appendix  B.  The 

aerodynamic  quantities  are  expressed  as  non-dimensional  coefficients 

and  derivatives  that  fit  the  data  package  (Appendix  A) .  The  resulting 

equations  of  motion  to  be  programmed  are  as  follows: 

2 

u  *  -g  sin  8  +  rv  -  qw  +  S-  (C  +  C  5e) 

-m  x  x5e 

+  ^  (T^  +  T2)  cos  otT  cos  3T  (8) 
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o 


o 


p 


I  I  -  I  2 
x  z  xz 


_5.  pV2Sb  [C  +  C  6a 

2  *  V 


+  C  6r+-(C,p  +  C,  r)  ] 

6r  p  r 

+  lz  [  (Tx  -  T2)  yT  sin  -  0^  +  T2)  zT  sin  6T 
I 

+  -M  PVzSb  (Cn  +  C^Ja  +  C^Sr 

+  W  (Cn  p  +  Cn  r)1  +  lKZ  [(T1  -  V  yT  co3  BT 
P  I 


-  (I,  +  y  ^  sin  6t]  +  Ix2  (wTi  -  wl2>  x„ 


+  (I  -  I  +1)1  pq  +  (I  I  -  i  2  -  I  2)  qr 
x  y  z  xz  y  z  z  xz  n 


+  I  I  q  co 
xz  r  1  r 


\  PV!S'  «„  +  c  **  +  fv  ■» 

oe  q 


-  (Tx  +  T2)  xt  sin  aT  +  (T1  +  T 2>  z T  cos 


+  (r2  -  p2)  +  CIZ  -  lx)  rp-rlr 


r  -  v 


y  pV2Sb  [C  +  C  6a  +  C  6r 
2  n  n6a  n6r 

+  2V  (Cn  p  +  Cn  r)]  +  (I1  '  T2>  yT  cos  61 
P  r 

“  (T1  +  T2)  Sln  BT  +  <\  -  V  X«T  +  Xxz 


(ID 


(12) 


(p  -  qr) 


+  (Ix  -  ly)  pq  +  q  ir  «r) 


The  Euler  relationships  are  developed  in  Appendix  B. 
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Control  System 

A  key  objective  of  this  study  is  the  comparative  evaluation  of 
the  response  of  the  aircraft  model  to  control  surface  deflections  with 
varying  degrees  and  types  of  chrust  augmentation.  For  this  reason, 
the  influences  of  the  F-lll  self-adaptive  flight  control  system  and 
the  dynamics  of  pilot  response  are  purposely  removed  from  the  control 
loop.  A  simple,  closed-loop,  feedback  control  law  based  on  yaw  (8), 
yaw  rate  (r),  and  roll  rate  (p)  is  used  to  command  aileron  and  rudder 
control  surface  deflections. 

«a  e  K!  P  (14) 


6  “  r  -  K,  (0  -  8  ,) 

r  2  3  command 


The  control  gains  are  purposely  kept  at  the  minimum  values  necessary  to 
maintain  wings-level  at  low  angles-of-attack.  The  low  gain  levels  will 
not  inhibit  departure  or  stall  in  the  high  angle-of-attack  regime  by 
inadvertent  over-control.  Maximum  control  deflection  angles  are 
specified  in  Appendix  A. 

Thrust  augmentation  authority  limits  (TVA'J)  shown  in  Figure  3  are 
used  to  schedule  both  the  engine  exhaust  deflection  system  and  the 
a'  xiliary  thruster  system.  Since  control  augmentation  is  necessary 
only  in  the  high  angle-of-attack  regime,  the  authority  is  scheduled  as 
a  function  of  angle-of-attack,  rudder  deflection  angle,  and  a  maximum 
specified  thrust  control  moment.  The  engine  exhaust  nozzles  are 
assumed  to  lie  in  the  x-y  plane  and  the  thrust  deflection  is  restricted 
to  produce  only  yawing  moments  and  longitudinal  and  lateral  forces.  To 
facilitate  comparison  between  the  two  systems,  the  maximum  moment 
produced  by  the  auxiliary  thrusters  is  identical  to  that  produced  by 
engine  exhaust  deflection. 
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III.  Simulation  Test  Plan  and  Evaluation 

The  simulation  and  evaluation  of  results  are  divided  into  two 
phases.  The  first  phase  is  run  solely  to  generate  dynamic  stability 
parameters  for  co'.^arative  evaluation  of  each  configuration.  The 
second  phase  consists  of  several  "flights"  of  each  configuration  at 
high  angles-of-attack  using  various  flight  control  parameters  to  pro¬ 
duce  maneuver  envelopes  based  on  criteria  to  be  defined. 


Phase  _I 

One  simulation  is  run  for  each  of  the  following  configurations: 

a.  Basic  configuration  with  no  thrust  augmentation 

b.  Engine  exhaust  deflection  augmentation 

(1)  50%  authority  limit  (6°  maximum  deflection) . 

(2)  100%  authority  limit  (12°  maximum  deflection) . 

c.  Auxiliary  thruster  augmentation 

(1)  50%  authority  limit 

(2)  100%  authority  limit 

To  provide  a  basis  for  comparison,  the  authority  limits  and  algorithm 
for  the  auxiliary  thrusters  are  designed  to  produce  moments  equal  to 
those  developed  by  the  engine  e.haust  deflectors  without  the  force 
effects. 


Evaluation  of  each  simulation  in  this  phase  will  be  based  on  the 
static  lateral  stability  parameters  (Cn^  and  C^)  and  the  dynamic 
directional  stability  parameter  (Cn^  ^  ) .  The  derivation  of  Cn^ 

can  be  found  in  Reference  19.  This  report  and  others  (Ref  5  and  28) 
show  good  correlation  between  divergence  characteristics  pr^dJcted  by 


^Bjdyn 


and  actual  aircraft  divergence  characteristics. 


The  method 
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for  determining  these  parameters  is  different  than  the  normal  method  of 
extracting  Cn$  and  Cj^  from  wind  tunnel  data  and  then  calculating  a 
^  .  Since  thrust  effects  do  not  appear  in  the  characteristic 

equation  used  for  the  derivation  of  Cng  ^  ,  terms  are  defined  to  com¬ 

bine  both  the  aerodynamic  and  thrust  effects  as  one  parameter.  The 
elevator  control  will  produce  a  gradually  increasing  angle-of-attack  up 
to  stall  or  departure.  Aileron  control  will  be  used  to  minimize  roll 
and  maintain  wings  level.  The  rudder  and  thrust  augmentation  devices 
will  be  driven  by  a  sinusoidal  commanded  sideslip  angle.  The  resulting 
motion  is  large  sinusoidal  oscillations  in  yaw  and  yaw  rate  with  a  mini¬ 
mum  of  variation  in  the  other  states. 

An  effective  Cn  and  is  calculated  at  each  time  interval  (0.1 
seconds)  in  the  simulation  with  the  thrust  effects  included. 


pV2Sb 


21 

_ z 

pV2Sb 


^thrust 


N 

thrust 


(16) 

(17) 


C£q  and  CnQ  are  interpolated  values  calculated  in  the  computer  simula¬ 
tion  based  on  the  data  (Appendix  A)  and  exclude  the  effects  on  the 
rolling  and  yawing  moments  produced  by  roll  and  yaw  rates. 

The  change  in  sideslip  (Aft)  is  calculated  for  the  same  time  incre¬ 
ment  and  the  static  lateral  stability  parameters  arc  estimated  by 


l 

ft.cff 

Cn 

ft.eff 


AC 


£.ef  f 


Aft 


AC 


neff 


Aft 


(18) 

(19) 
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and  the  dynamic  directional  stability  parameter  is  estimated  by 

I 

C  -  C  cos  a  -  “  C.  sin  a  (20) 

°3,dyn  n3,eff  x 

The  values  of  effective  Cn^,  »  and  Cn^  are  determined  and 
plotted  as  a  function  of  angle-of-attack.  The  results  are  an  indica¬ 
tion  of  the  combined  aerodynamic  and  thrust  effects  on  the  stability  of 
the  system. 

Phase  II 

For  each  simulation  in  this  phase,  the  aircraft  configuration  is 
flown  to  departure  or  stall.  To  provide  variation  in  departure  atti¬ 
tudes,  sideslip  command  (3C)  in  the  rudder  control  (Equation  15)  will 
be  set  at  a  different  value  for  each  simulation.  Several  simulation 
runs  are  made  for  each  of  the  following  configurations: 

a.  Basic  configuration  with  no  thrust  augmentation 

b.  Engine  exhaust  deflection  augmentation 

(1)  50%  authority  limit 

(2)  100%  authority  limit 

c.  Auxiliary  thruster  augmentation 

(1)  50%  authority  limit 

(2)  100%  authority  limit 

d.  Engine  exhaust  control  with  rudder  fixed  at  neutral 

. (1)  50%  authority  limit 

(2)  100%  authority  limit 

In  the  final  sets  of  simulations  (Part  d),  the  rudder  will  be 
fixed  at  ?.ero  deflection  and  the  deflected  engine  thrust  will  substi¬ 
tute  for  the  rudder  as  a  lateral  control  device.  The  authority  for 
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each  of  ti  -se  configurations  will  use  the  limits  indicated  for  all 
angles-of-attack.  The  thrust  control  will  be  maintained  at  full 
authority  limits  for  all  angles-of-attack.  The  purpose  of  the  simula¬ 
tions  in  Part  d  is  to  provide  so  me  comparison  of  the  control  effective¬ 
ness  of  rudder  alone  versus  thrust  alone. 

To  evaluate  and  compare  the  resulting  time  histories  of  each 
simulation,  divergence  boundaries  are  established  based  on  two 
different  criteria.  The  first  compares  the  sign  of  the  first  and 
second  time  derivatives  of  beta  (S) .  If  each  is  of  the  same  sign  (both 
positive  or  both  negative),  the  time  and  aircraft  attitude  is  recorded 
and  a  departure  condition  is  defined.  The  series  of  departure  condi¬ 
tions  recorded  for  each  configuration  are  plotted  on  a  graph  with 
angle-of -sideslip  (8)  as  the  ordinate,  and  angle-of-attack  (a)  as  the 
abscissa  and  labeled  the  "Beta  Condition  Departure  Envelope".  The 
second  criterion  of  departure  is  the  completion  of  20°  of  spin;  when 
this  occurs,  the  aircraft  attitude  and  time  at  the  20°  point  are  re¬ 
corded.  The  results  are  plotted  in  a-S  coordinates  and  are  labeled 
"Spin  Condition  Departure  Envelope". 

It  is  important  to  realize  that  these  envelopes  are  not  the 
"maneuvering  envelopes"  defined  in  the  flight  manuals.  Both  are  simply 
diagnostic  tools  for  comparing  the  effectiveness  of  each  thrust  aug¬ 
mented  configuration. 
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IV.  Results 


Phase  I 


The  first  simulation  of  this  phase  shown  in  Figure  4  was  run  with¬ 
out  thrust  ajgmentation.  The  resulting  Cn^  was  compared  with 
previous  results  (Ref  19  and  2d)  and  provided  baseline  data  against 
which  to  evaluate  the  effectiveness  of  the  thrust  augmented  control. 
Correlation  was  limited  to  angles -of-attack  less  than  35°  which  was 
the  maximum  angle-of-attack  in  the  previous  calculations.  The  present 
method  gave  surprisingly  good  agreement  and  established  confidence  in 
the  method  used  in  the  computer  program  for  calculating  an  effective 
directly  from  the  simulation.  The  remainder  of  the  Phase  I 
simulation  results  are  shown  in  Figures  5  through  8.  The  apparent 
flattening  of  the  curves  in  Figure  5  is  due  to  a  scale  change. 

The  most  significant  effect  of  thrust  augmentation  was  the  exten¬ 
sion  of  the  angle-of-attack  boundaries  of  departure.  The  fact  that 
simulations  did  not  depart  when  the  calculated  Cn^  ^  became  negative 
was  probably  due  to  the  simplified  control  system  used.  As  shown  in 
Table  1,  both  the  sign  change  cl  C;1^  ^  and  the  maximum  attained  angle- 
of-attack  are  improved  with  the  thrust  augmented  configurations. 

Since  the  authority  limits  for  thrust  control  were  based  on  the 
rudder  effectiveness  (Cn,$r)>  it  was  noted  that  the  lower  limit  selected 
occurred  right  at  the  Cng  crossover.  At  this  point,  the  authority 
limits  were  altered  to  allow  the  thrust  augmentation  to  begin  at  an 
anglc-of-attacl.  of  20°  and  reach  full  authority  at  an  angle-of-attack 
of  40°.  The  simulations  were  run  again  and  the  results  (Appendix  C) 
indicated  decreases  in  the  angle-of-attack  where  the  sign  of  the 
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FIGURE  4.  Lateral  Stability  Characteristics 
(No  Augmentation) 


FIGURE  5.  Lateral  Stability  Characteristics 

(Exhaust  Deflection  with  5C%  Authority  Limit) 


FIGURE  6,  Lateral  Stability  Characteristics 

(Exhaust  Deflection  with  100%Authority  Limit) 
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effective  C™  ,  became  negative.  This  indicated  that  rather  than  an 
np, dyn 

increase  as  expected,  the  lateral  stability  of  the  thrust  augmented 
configuration  had  decreased.  For  verification,  the  simulations  were 
again  run  with  complete  time  histories  and  the  departure  angle-of- 
attack  for  each  thrust  augmented  configuration  with  the  20°-40° 
authority  limits  was  approximately  5®  lower  than  that  of  the  30°-50® 
authority  limits.  Since  the  original  authority  limits  with  initial 
thrust  deflection  at  30°  angle-of-attack  and  maximum  thrust  deflection 
at  50°  angle-of-attack  had  produced  satisfactory  results,  it  was 
decided  to  proceed  to  Phase  II  with  the  authority  limits  unchanged. 


TABLE  1 


Summary  of  Phase  I_  Simulations 


Configuration 


Maximum 

Deflection 

Angle(deg) 


An gle-of -Attach 
^B.dyn 

Sign  Changa(deg) 


Angle-of-Attack 
at  Departure 

(deg)  _ 


No  Augmentation  0 
Exhaust  Deflection  6 
Exhaust  Deflection  12 
Auxiliary  Thrusters  *  6 
Auxiliary  Thrusters  *12 


30.4 

50.0 

31.7 

55.6 

30.7 

61.5 

30.7 

55.6 

30.8 

61.5 

*  Equivalent  in  moment  to  exhaust  deflection  of  maximum 
angle  indicated. 


Phase  II 

The  intent  of  this  phase  was  to  provide  an  alternate  method  of 
evaluating  the  effects  of  thrust  augmentation,  independent  of  the 
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method  used  in  Phase  I.  The  rationale  for  this  method  arises  from  the 
apparent  relationship  which  exists  between  angle-of-attack  (a)  and 
sideslip  (3)  at  the  time  of  departure.  The  idea  was  to  select  3  and 
Increase  a  until  departure  occurs  over  a  wide  range  of  sideslip  in 
order  to  form  an  envelope  which  would  describe  the  non-linear  nature  of 
their  relationship  at  departure.  Unfortunately,  the  piocedure  proved 
to  be  time  consuming  and  costly  in  relation  to  the  value  of  the  results 
obtained. 

The  key  to  defining  the  envelope  was  the  development  of  proper 
criteria  to  accurately  determine  the  departure  condition.  Using  the 
first  criterion,  the  first  and  second  time  derivatives  of  sideslip,  it 
was  found  that  departure  occurred  in  the  first  seconds  of  simulation 
while  examples  of  the  full  time  histories  shown  in  Appendix  C  show  the 
departure  occurring  at  times  of  20  to  30  seconds.  Several  variations 
in  the  method  of  calculating  the  rates  and  applying  the  test  were  tried 
in  an  effort  to  improve  the  criterion.  For  example,  minimum  magnitudes 
of  3  and  3  or  minimum  elapsed  times  were  imposed  as  conditions  required 
before  the  test  could  be  applied.  None  were  successful;  the  computer 
program  indicated  false  departure  immediately  after  the  imposed  limits 
were  exceeded. 

Similar  difficulties  were  encountered  using  the  developed  spin 
condition.  Since  the  angles-of-attack  and  sideslip  oscillate  violently 
and  rates  are  unpredictable  once  the  spin  has  developed,  it  was  diffi¬ 
cult  to  predict  how  far  to  backtrack  from  the  test  point  to  the  actual 
departure  point. 

Examination  of  the  time  histories  did  not  indicate  why  these 
methods  failed.  The  correct  criterion  for  defining  departure  will  have 
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to  be  based  on  several  of  the  simulated  states  due  to  the  non-linear 
nature  of  the  problem. 

Overall  Results 

Sample  time  histories  of  simulation  runs  for  some  of  the  configur¬ 
ations  are  shown  in  Appendix  C.  With  the  control  law  used,  the  model 
in  each  case  increases  in  angle-of-attack  until  departure  occurs. 
Examination  of  full  time  histories  i:1  licated  that  with  thrust  augmenta¬ 
tion,  the  maximum  angle-of-attack  reached  before  departure  is  increased 
by  5°  to  10°.  There  is  little  apparer :  increase  in  the  g  loadings;  all 
are  well  within  pilot  tolerances. 

Surprisingly,  the  use  of  lateral  thrust  control  with  the  rudder 
held  fixed  provided  controllability  nearly  equal  to  that  of  the  rudder 
alone.  Obviously,  the  vertical  tail  and  fixed  rudder  have  a  signifi¬ 
cant  stabilizing  effect.  This  simulation  does  not  infer  the  use  of 
thrust  in  place  of  the  rudder,  but  suggests  that  thrust  deflection 
could  provide  emergency  control  device  if  the  use  of  the  rudder  is  lost, 
a  condition  quite  possible  during  or  after  combat  engagement- 

Time  histories  for  some  of  the  simulations  used  in  Phase  I  ave 

also  shown  in  Appendix  C.  The  magnitude  of  the  sinusoidal  sideslip 

command  ($c)  is  unrealistic  for  a  real  aircraft,  but  provides  a  suitable 

and  simple  method  for  extracting  C  ,  Cp  and  C 

n8  6  |3,dyn 
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V.  Conclusions  and  Recommendations 


Conclusions 

The  results  of  the  Phase  I  simulation  indicate  that  both  the 
engine  exhaust  deflection  and  the  auxiliary  thrusters  are  effective 
lateral  control  augmentation  devices.  The  maximum  angles -of -at tack 
attained  prior  to  departure  using  either  of  the  augmented  configura¬ 
tions  were  up  to  11.5°  greater  than  those  of  the  conventional  model. 

The  calculation  of  an  effective  Cn^  compared  well  with  values 
from  a  previous  study  (Ref  19)  for  the  same  aircraft  configuration. 
Although  it  could  not  be  used  to  predict  the  exact  departure  point  in 

the  model  used,  the  effective  C_  ,  was  a  useful  tool  for  selecting 

n$,dyTi  6 

the  thrust  control  law  authority  limits. 

The  presence  of  side  forces  and  th:  slight  reduction  in  longi¬ 
tudinal  thrust  present  with  the  use  of  exhaust  deflection  had  little 
effect  on  the  response  of  the  model.  The  predominate  control  factor 
was  the  moment  generated  by  either  augmented  configuration. 

Thrust  deflection  alone  was  shown  to  be  an  effect  device  for  main¬ 
taining  control  in  the  event  of  rudder  control  loss. 

Recommendations 

The  following  recommendations  are  made  regarding  the  use  of  thrust 
control  augmentation  in  high  performance  aircraft. 

1.  An  appropriate  follow-up  to  this  study  would  be  to  investigate 
the  hardware  aspect  in  detail  to  determine: 

a.  expense  of  implementation 

b.  time  responses  of  various  control  hardware  (considerable 
work  has  been  accomplished  in  the  areas  of  thrust  control  on  V/STOL, 
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aircraft  and  RCS  rockets  on  re-entry  vehicles). 

2.  The  use  of  augmented  control  by  thrust  vectoring  should  be 
applied  to  other  multiple  and  single  engine  aircraft  with  histories  of 
lateral-directional  stability  problems  at  high  angles-of-attack. 

3.  A  suitable  criterion  for  determining  the.  departure  point  in 

the  six-degree-of-freedom,  non-linear  simulation  needs  to  be  established. 
This  would  allow  the  calculation  of  an  a-3  maneuver  envelope  from  the 


computer  simulation. 
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Appendix  A 
Aircraft  Model  Data 
(Ref  20  and  23) 


A  three-view  of  the  F-111A  modeled  for  this  evaluation  is  shown 
in  Figure  A-l.  The  non-linear  aerodynamic  coefficients  and  derivatives 
are  presented  in  Table  A-l.  The  following  are  the  aircraft  model 
general  parameters. 

Overall: 

Length . 72.13  ft 

Height. . 17.12  ft 

Weight . 50,000  lbs 

Wings : 

Span . 63.0  ft 

Area . 525  ft2 

Sweep . 26c 

Mear  Aerodynamic  Chord . 9.04  ft 

Aspect  Ratio . 6.97 

Inertial  Terms: 


Ix . 50,000  ft2  -  slug 

Iy . 315,200  "  " 

I  . 351,500  "  " 

z 

i  .  o  "  " 

xz 

I  .  0  "  " 

r 
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Maximum  Control  Surface  Deflections: 
<Se  10,  -  25  deg 

6  ±15  deg 

±30  deg 
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Appendix  b 
Equations  of  Motion 
(Ref  14,  15,  16,  and  25) 


Iv  ' 


Assuming 

1)  negligible  earth  rotation, 

2)  negligible  changes  in  moments  of  inertia, 

3)  x-z  is  a  plane  of  symmetry, 

4)  h^  is  the  only  significant  rotor  term,  and 

5)  negligible  changes  in  the  rotor  term, 

the  force  and  moment  equations  in  the  body-fixed  reference  frame 
(Euler's  equations)  are 

X  -  mg  sin  0  =  m[u  +  qw  -  rv]  (B-l) 

Y  +  mg  cos  0  sin  <J>  ■  m[v  +  ru  -  pw]  (B-2) 

Z  +  mg  cos  6  cos  <{>  =  m[w  +  pv  -  qu]  (B-3) 

i*  ■  -  \z  (r  +  pq)  -  <xy  -  xz>  qr  (b-4) 

M  *  V  "  (r2  "  p2)  "  (IZ  -  y  *p  +  r  IrWr  (B~5) 

N  -  Izr  -  Ixz  (p  -  qr)  -  (Ix  -  Iy)  pq  -  ql^  (B-6) 

Rearranging  the  equations  for  use  in  the  computer  program, 

u  “  -  g  sin  0  -  qw  +  rv  (B-la) 

m 


v  .  _  4.  g  cos  0  sin  <J>  -  i.  u  +  pw  (B-2a) 

m 


w  ■  —  +  g  cos  0  cos  $  -  pv  +  qu  (B-3a) 

m 


( 
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Simulation  Results 
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FIGURE  C-l.  Lateral  Stability  Characteristics 
(Exhaust  Deflect ion/TVAU  a  50%/Lower  Authority  Ranges) 


FIGURE  C-2.  Lateral  Stability  Characteristics 
(Exhaust  Def lection/TVAU  =  100%/Lower  Authority  Ranges) 


RirHR  -  DEC 


RLrHR  -  OEG 


RUN  NUMBER  -43M5 

FIGURE  C-3.  Lateral  Stability  Characteristics 
(Auxiliary  Thrusters/TVAU  =  502/Lower  Authority  Ranges) 
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FIGURE  C-7.  Simulation  Time  Histories 
(Rudder  Fixed/Thrusc  Control/TVAU  =100%) 


03S/C30  -  t  03S/030 


reuse  t  tv  -  rr/sec  aileron  otsp  -  ceo  Rt;30ER  otsp  -  oeo  elevator 


CftOM  BftWOC  -  rt  »10*  THRUSt  -  UB3  <10*  7  RCCCLERRrtON  -  O  T  ACCELERATION  -  0  X  ACCELERATION 

0<«tn  »«,«  ».»»  m  •«  _-i.OT  ~».i»  -o.ig  _-<.eg  «*.ia  o.«r  _p.70  o.« 


oeo/sec  e£TA-oor  -  oeo/sec 


f IGURK  09  (Cout.).  Sinulat ion  Ti:r.c  His"orics 
(Enfi;.ne  Deflection/TV, MJ  =1  )0%/Sinuooidal  8  ) 


CROSS  RANOE  -  FT  »tO*  fWT'JST  -  LB3  «10*  7  ACCELERATION  -  0  Y  ACCELERATION  -  0  X  ACCELERATION 

-«9,«9  9,50  J.00  J9.B9  *0.00  _*3 .  33  *1.39  l  .00  -J-3.S9  *0.10  0.30  J.30  0.48 


GAM/ AE/ 7 3-14 


Appendix  D 


Computer  Program 


70 


oo  o  o  o  o  nnnnnno  ooooononoooortortoono 


GAM/AE/ 73-14 


PROGRAM  RESFCNSdNF-UT  ,  CUTPL'T>T  AFE5=I  N°UT  »TAPE6=CUTPUTtPL0T  ,TAFE7) 


***«********>r*************  **********************************  ****** 
***««**««**««*«*«*****«*.*««**«****«*****************«************ 


THIS  PROGRAM  IS  A  SIX  DEGREE  OF  FREEDOM  SOLUTION  USING 
THE  GENERAL  NON-LINEAR,  BODY-AXIS  SYSTEM  EQUATIONS  CF 
MOTION  ANC  EULER  ANGLES... 

OEAS  M .  HARLEY,  2LT,  US  AF  AFIT/EN/  GAM-73 

**************  **************  **********  «**  *******************  ****** 
**************«***********;  *********************************  ****** 


**♦  DIMENSION  STATE  VARIABLES 
DIMENSION  Y  M2),  YP  (12) 

***  COMMON  FILES  NAMED  FOR  SUBROUTINES  AND  MAIN  PROGRAM 


1)  MAIN  PROGRAM  ANC  SUBROUTINE  GENERAL  PARAMETERS 

CCHMCN/CmAIN/EFCMO.AV,  AVC5  ,AVOT,JA  ,AO,  PI,  TOOT,  TRIMCFT,  NEAR  , 
1THRUS1,TH3L'SE,A1,A2,  A  3, 01,  62,3  T  ,T  1  ,T  2,  F  EG ,  S IGHC  ,  THE  C ,  ANGLE  , 

2GA, GP,GC*Cc*CA, OR, TRTl,TrTZ,AL,6E,ALPHAT, START, CNEE  ,CLEE,VCFT, 
6SFE,CFE, VR.CTCR, ONER, CHOP, THRU S ,T» IM , BE T AT  ,  WTl , WT2 , XWT ,C ON T , ZMT  , 
7AL0T, SECT, ALCFT,3E OPT, TCPT ,XT,  YT,ZT,  IX, IY,IZ,IXZjIR,WR,MA$S»G,S, 
8B,C»CX,CX0E , CZ,C70£tCM , CMO E ,CMQ ,RV R2 ,ROF , S T HE , C THE , TVGAIN 

2)  AIRCRAFT  AERCCYNAMIC  COEFFICIENTS  IN  NCN-C3MENSIO NAL  FORM 

CCMMCN/CCCEF/FCY (2 1, 9 ) ,ECN (21,  9 ) ,ECL (21 , °> ,ECM ( 2 1 , 9 ) , ECX (2 1 ) , 
1ECXDF  (21)  *  Fr  2 C E (  21 ) ,  ECMCH  21) ,  ECHO  (2  1)  »  ECYGR  (2 1 )  ,EC NOP  (2 1)  , ECLC R  ( 2 
21)  *ECYCA(21)  ,ECNOA  (21)  ,ECLCA(2l)  ,ECYP(21)  ,ECNP  (21)  ,  ECLF(21  )  » £CYR  (2 
311»ECN’R(21)»ECLR(21)  ,ECZ(21) 

3)  PARAMETERS  USFC  BY  PLOTTING  SUBROUTINE 

CCHMCNV/ALFHA  (40  2)  ,eETA(40  2),VEL(4  02  ),P(4C2),LCCF,0(A0  2), 

1R  ,  .02) , ALT (A  (E), THETA  (AG2)  ,PHI ( h92)  ,  PSI (402 ) ,YC  (7) , TIME( AO  2)  , 

2D  1ST (402) ,TRAV(402) >  AX  (4 02 ), 4* (40?) , AZ ( 402) , DELE (40 2) ,CELA  (402)  , 
3CELR(4P2) ,  TG\  (40  2)  ,THRST(4C2) ,  CNE(40  2)  ,CLE  (4C2)  ,CNGEFF  (4  02  )  , 
4LOP,CALF(4C2),ALOOT(  402)  ,RFOOT  (  402) 

C 

C 

EXTERNAL  GYRATES 
REAL  IX,  IY  ,IZ,1XZ,  IR^ASS 
C 

c 
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*•*  INPUT  TPE  TRIM  OPTICS  (TRIMOPT),  CONTROL  GAINS,  AIRCRAFT 
GEOMETRIC  ANO  CONTROL  PARAMETERS,  AIRCRAFT  AEROCYNAMIC 
COEFFICIENTS ,  ANC  INITIAL  CONDITIONS. 

N  AMELIST/C  AT  AS/TRI  SOFT,  GE,  CA,DR,Y0 ,S IGH 0  ,THEC , FEC ,L COP 
NAMELIST/GEC  N/IX,IY,I7 ,1X2 , IR,  WR,M AS S ,S , C , E , THFUSl , TMRL32, XT , YT  , ZT 
NAMELIST /CCEF/ EC Y, CCN , ECL* ECM,  ECX, EC  7 , E CXCE , FCZCE , ECMOE, EC MC,ECY0R 
1 , ECNOR , ECL  CF ,ECYDA ,ECLCA,ECNOA ,£CYP, ECNF ,ECLF, ECYR, ECNR, EC LR 
READ (5  ,CAT  AS  ) 

READ(S,C-ECn 

REAC(5,CCEF) 


DEFINE  PARAMETERS  USED  AS  CONSTANT  COEFFICIENTS  IN  THE  STATE 
EOLATIONS  AND  AS  ANGLE/RADIAN  RATIOS. 


ALPHAT^PETAT  =  7  =  0 . 
G-32. 1740 
A«OT=.  1*AV 
TCPI=2.*FI 
A2=S*P/2. 
B1=IX7MIY-IY  +  I?) 
T 1=THPUS 1 
TRT1=T1/THFUS 


3  NEA  R=0 
3  AO=l./(2*PI) 

3  A  VO  5=2.  *  AVOT 
3  WTl=WT2  =  XWT=ZMT=0 


S  PI  =  3  •  141592  653  E 
t  AV=180./PI 
t  ANGLE-PI/2. 

S  Ai=S/(2.*MftSS) 


3  A3=S*C/2.  ?  READER=TRIM  =  CHCP=S  TAPT  =  C. 

$  E2=  I7*1Y-IZ**2-IX7**2  ?  B3=IX*IZ-  IX  ?**  2 
3  T  2=  THRU  S  2  S  7MRUS  =  TH-T2 

3  TRT  2=T2/THRUS  3  XKT  =  5C. 


***  DEFINE  INITIAL  CONCITIONS. 

0060 J=1 ,7 
60  Y(J)=YC(J) 

Y(0>=THEO?Y(S)=SIGHOJY(1G)  =FE0?Y(11)  =Y(12)=0. 


***  SELECT  CESIFEO  TR Iv  OFTIONS  (TRIHOPT)  USING  THE  FOLLOWING 
COCE  l 

0.  -  NO  TRIM 

1.  =  STRAIGHT  ANO  LEVEL 

1.  -  STEADY  CL  INC  (WITH  OROPEP  INITIAL  CONCITIONS) 

1,  -  STEAOY  CESCENT  (WITH  PROPER  INITIAL  CCNOITIONS) 

2.  =  COORDINATED  TURN 


IFITPIMCFT.EC.C.)  GO  TC  9 
IFCTRIKCFT . E  C . 2. )  GO  TC  10 


***  TP IM  CCNC1TION  FOP  STRAIGHT  AND  LEVEL  FLIGHT,  STEADY  CLIME,  CR 
STEACY  CESCENT 

CALL  TPIM1  (Y  ,YF) 

IF(CHCF.EC.i  .)  GO  TO  6 
GC  TO  9 
10  CONTINUE 
C 

/2 
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***  TRIM  C0NC1TI0N  FOR  TURNING  FLIGHT. 


CALL  7RIH2  (Y  ,YF) 

IF (CHCF. EC . 1 .) GO  TO  6 
9  CONTINUE 

***  PRINT  CL T  FIXEO  PARAMETERS  AND  INITIAL  TRIMMED  CONDITIONS 

CALL  PRINT  1  (  Y ) 

Y (2) =1. 

***  THE  SOLUTION  IS  NOW  WORKED  USING  THE  RKGXY?  INTEGRATION 
ROUTINE  ANC  THE  DESIRE  C  ATPCRAFT  CONTROL  SUBROUTINE. 

0077J=1,LCCF 
CALL  CCNTRL  (V,YP,T) 

IF(CCNT.EC.O.)C-0  TO  5 
CALL  VECTOR ( Y ,YP) 

5  CONTINUE 
00  76  JFK=1, 2G 

CALL  F  KG  XY  Z(  T,Y,  YP  ,  12,  .  0  G5  ,.0D  COO  0  5,  GYRATES) 

IFlN'EAR.EG.l  )  GO  TO  87 
76  CONTINUE 

OESIREC  PARAMETERS  ARE  STORED  FOR  THE  PLOTTING  SUBROUTINE 
GA=J 

CALL  STORE  (Y,YF) 

77  TIME ( J ) =T 
GO  TO  68 

87  LCOc=J-l 

88  CONTINUE 

**»  SUERCUTINE  CUTPLOT  IS  CALLED  AND  RESULTS  ARE  FRINTEC 
AS  A  FUNCTION  OF  TIME. 

CALL  SPEVAL(Y,Y°) 

CALL  CLTFLC T 

6  CONTINUE 
STOP 
END 
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SUBROUTINE  ORATES  (T  ,Y  ♦  YP) 


***  THIS  SUGR  CLTINE  CONTAINS  THE  STATE  EQUATIONS  USEC  BOTH  IN  SETTI 
THE  TBIM*EC  CONDITION  AND  IN  THE  PROGRAM  SOLUTION.  DATA  IS 
LINE  AR 12  EC  WITH  RESPECT  TO  THE  PARTICULAR  AIRCRAFT  ATTITLCE  EA 
TIME  THIS  SUBROUTINE  iS  CALLED.  IT  IS  CALLED  BY  THE  MAIN 
FROGRAH  A  NO  AS  AN  EXTERNAL  BY  SUBROUTINE  RKGXYZ • 


OIHENSICN  Y ( 12) ,  YP (1 2 ) 

COMHCN/CHAIN'/BECMD  ,AV,  AVC5  , AVO T , JA , A  0, P I , TCPI ,  TRIM  OFT,  NEAR  ; 
1TWRUS1,THRUS  2,fl, A2, A3  ,E1, E2,B3  ,T1 ,T 2 , F E C , S IGHC , T HE C , ANGLE , 
2GA»Ge,GC» CEtCA*3R. TRT1 ,TRT 2, Al , PE, ALPHA T, START, CNEE »CLEE,VCPT  , 
6SFE,CFE,VR  ,  C  YCR,  CN  CR  ,  CHCF,  THR'JS  ,TR  IM  ,  PET  AT  ,  HT  1 ,  NT  2 ,  X  WT  ,CON  T  ,  Z  »  T  , 
7ALDT,GECT,ALCFT, 3ETPT ,TCFT ,XT,  YT,ZT,  IX , I Y , IZ , I XZ , IR , HP ,MAS  S  ,  G  , S  , 
fiB,C,CX,CXCE,CZ,CZDE,CM,CH.DE,CM0,RVR2  ,  ROH  ,  STHE  ,  C  THE  ,  T  VGAI N 

COMHCN/CCOEFyECY (2 1, 9) ,ECN  .21, 9), E CL (21, B) ,ECM(2l,9),ECX(2 1) , 
1ECX0E  (21)  ,  EC2CE121 )  ,ECHCE(  21),  ECHO  (21)  ,ECYCR(21)  ,ECNCR  (21)  ,ECICF  (2 
21) , ECYCA (21)  f ECMOA  (21)  ,  ECL DA (2 1 ) ,ECY P ( 2 1 ) , ECNP ( 21 ) , ECLF ( 21 ) , E C Y R ( 2 
31),ECNR(21)  ,ECLR(21)  ,  EC  7  <2  1  > 

■REAL  TX,IY,IZ,IXZ,  IR.MASS 


VARIABLE  COEFFICIENTS  ARE  OEFINEO  FOR  USE  IN  THE  STATE  EQUATION 


VRs=SQRT(Y(  1)  **2h  Y (  2)  *  *  2* Y(  2)  **2) 

RCH=.  00237  8*  (Jl,-,  00  0  0 06  85  * Y<7)  >**4.  256) 
FVR2=FCHM  VR«*c) 

VCVR  =  Y  (2  )  /  V  R 
CTHE=CCS  (Y  (9  )) 

IF(PBS(CTHE)  ,LT,  .0001)  NEAPrl 
IF(NEAR.EO.l)  RETURN 
IF(A0S(Y(1))  .LT.1.E-2C)  GO  TO  53 
AL=ATAN2(Y  (2  )  ,Y< 1)  ) 

GC  TO  54 

53  AL  =  SIGN  (ANCLE, Y( 3) ) 

54  CONTINUE 
BE=ASTN  ( VCVR  ) 

SFE  =  SIN(Y  (1  C  ) ) 

STHE=SIN (Y  (S  )) 

CFE=CCS  (Y  ( 1  C  ) ) 

SPSI  =  SIN  (Y  <<5 ) ) 

CFSI  =  CCS(Y  <<=)) 

SET=SIN(BFTAT) 

C3T =CCS ( RE  T  A Y ) 

S A T  =  S I N  (ALFHAT) 

CAT=COS(ALFH£T) 
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***  AERCOYNANI C  COEFFICIENTS  ARE  EXTRAPOLATED  AS  A  FUNCTION 
OF  THE  AIRCRAFT  ATTITUCF  (ALPHA  AND  PETA). 

JA=AVC5*AL*3 . 

JE=AVCT*PE*S . 

IF((Al*AV)  .CE.Sa.) JA  =  2  0 
IF((AV*AL)  .IE. -10.  )JA  =  1 
IF((AV*0£) «L  E • -4  0  • )  Jn  =  1 
IF ( ( A  V* EE )  .C-f  .40.)  J9=6 
IALI=S*(JA-2) 

ALI=I ALI 

0109=  (AV*AL-ALI) /5 * 

IPEI  =  10*  (JE-E) 

8EI-IGEI 

OELOT=(AV*EE-FEI)  /10« 

CX-(ECX(JAM  )-ECX  (  JA  )  )  ’CLOE4-ECX  (JA) 

C2=(EC2(JA4l  )-EC7(JA)  )  *CL0  E+EC  Z  (JA  ) 

CXDE=(ECXCE  UA+1) -ECXCE  <  JA  ) )  *DL  C5*-  EC  XOE  (JA) 

CZDE=(ECZCE  (Jf  +1)  -  EC  ZEE  (JA  ) )  *0  L  C5*  EC  ZDE  (JA) 

CMOE=(£CHCE  (JA+D-ECHCE  (JA  ) )  *nL  C5«-  EC  *DE  (JA) 

CPO  =(£CMQ  (JA  +  D-ECHC  ( JA  ) )  *0  L  CP*  EC  HO  (JA) 

CYl)R  =  (ECYCR  (wA*l)-fCYCP  (JA)  )*3LC5+EC  YOR  (JA) 

CNOR=(ECNCK  (*.A«t)  -  EC  NCR  ( J'  M  *0105*  EC  NOR  (JA) 

CLOR=(ECLOF  (  JA+1) -ECLCR  (JA  ) )  *0L  CO* EC  LOR  (  JA) 

CYOA=(ECYCA (wA  +  1) - ECYC A  (JA ) ) *0 L CP* EC  YD  A  (JA) 
CNDA=(ECNCA(wA+l)-CCNCA  (JA  )  )  *DL  C5+  EC  NO  A  (JA) 

CLOA=(ECLCA(  VA  +  1)  -ErLOA(JA  )  )  *0  L 0F+  EC  LD  A  ( JA  ) 

CYP  =<ECYF  (JAU)-ECYF  ( JA  )  )  *0  L  C5*  EC  YP  (JA) 

CNP  =(ECNF  <  J  A  41 )  -  EC  N  F  (JA  )  )#0LC5+EC  NP  (JA) 

CLP  =(ECLP  <  *.  A  «1 )  -  ECLF  ( JA  )  )  *1 L  05*  EC  LP  (JA) 

CYR  =(ECYR  (.A+D-ECYR  ( JA  ) )  *DL  C5+  EC  YR  (JA) 

CNR  ”(ECNR  (JA+D-ECNR  ( JA  )  )*D  L  C5*  EC  NR  (JA) 

CLR  =(ECLR  (JAU)-ECLR  ( JA  )  )  *0L  05*  EC  LR  (JA) 

CYAMECY  (JA  +  J,J3)-ECY  (JA  ,JE))*0LCP+ECY  ( JA,JP) 

CLA=(ECL  (JA*  J  ,J-3)  -ECL  (JA  ,JE))*0L05+ECL  (  JA,  JP) 

CNA=(  ECN  ( JA*  1  ,J0)  -ECN  ( JA,JP)  )*  DLOO  +  E  C‘l(  JA,  JR) 

CHA=(ECH(JA4  J,v3)-FCH(  jA,je))*CLC5>ECY(jA,„G) 

CYAF=  (ECY(JA4  1  .J34-1)  -  ECY  (J  A,J3  4l) )  ’0LO5  +  ECY  (JA  ,  JG  +  1) 
CLAP=(ECL(  JA  4l,JP4l)  -  ECL  (J  A  ,  J3  4  1) )  #0  L05  +  ECL  ( JA  ,  JR4 1 ) 

CNAP=  (ECN  ( JA  4  1,334  1)  -  ECN(JA,J8 4-1) )  ♦QL05+ECN  (JA  ,J341 ) 

C  M  A  P  =  (ECH(  JA  4l,JR4l)  -ECMJ  A,  J34j) )  *DL05  4ECH  (JA,  je4l) 
CY=(CYAP-CYA ) *CRLO  T  +  CYA 
CLs(CLAP-CLA  )*C3L0T4CLA 
CN=(CNAP-CNA  )»C3L0T+CN A 
CH=(CHAF-CHA  )  ’C3L0T4CHA 

Q1=CX4CXDE’CE 

C2=CY-»CYCA*C  A4CY0R»0R4  ,5*B*(CYP*Y(4)  +CYR*Y (61 > / VF 
03=CZ4C?nE*CE 

04=CL4CLCA*C  A4rL0°*CR4  .6*°* (CL  F»Y (4)  4CLR*Y  (6)) /VR 
GF=CK4CHCE*Cf  4.5*n*Cf'CrY  (5  )/7D 

Q6=CN 4CNC  A  *0  A  ♦  CNOR  ’OR  ♦  «  5  *R  *  (CN P*Y (  4)  4CNP’  Y  (6)  )  /  VR 

c 

C 
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444  THE  STATE  ECUATIONS  FOLLOW  WITH  THE  STATE  VARIABLES  DEFINEC  AS« 

Y  ( 1 )  =  L  Y  (  2)  =  V  Y (3 ) =  W  ECCY  AXIS  SYSTEM 

Y  C  4 )  =F  Y  {  5)  =C  Y  (6  )  =R  BCCY  AXIS  SYSTEM 

Y(C)  =  THETA  Y  (  9)  =  R  S I  Y(10)=PHI  EULER  ANCLES 

Y<11)-XE  Y  ( 12 )  =  YE  YC7)=ZE  EARTH  AXIS  SYSTEM 

YP<1) =-G4STHE4Y<6>  4Y  <2)-Y<  5>4Y ( 3) ♦ A1 4RV R24G If ( T1 4T2 ) 4CAT4C eT/HASS 
i* (HT 14-WT2) /MASS 

YF<2)=G4CTFE4SFE*r  (4 )  4  Y  < 3)  -Y  (6  )  4Y  ( 1)  MT1  +  T2)*SET/MASS 
liFVR24Al4G2 

YF(3)=G4CTHE4CFE*Y<5)  *  Y  ( 1)  -Y  (4 )  4Y  (  2)  -  ( T  1  +  7  2)  4S  AT/MASS 
1HRVR24 A  14C  2 

YP(4)=(I24  (R  VR  24  A2  404  4  (T 1-  T2)  *  S  AT*  YT  -  (T  1  4T 2 )  4SET4ZT  )  4-1  XZ4<  FVR24A24 
1C6-(T14T2)  *XT*S3T4  'T  1  -  T  2  )  *  YT4C  PT)  4  Bi  4Y  (  4 )  *  Y  (5 )  +  B24  Y  (5 )  4Y  (6  )  ♦  Y  (5  )  *  I 
2R*WR»IXZ)/E2  4  (NT  1-HT2)4XWT4IXZ/B3 

YP(5)  =  (RVR2* A3405-  (T 14T2 )4 XT4S A T4 { Tl +T2 ) 4ZT4CAT 4 IXZ4  <Y (6 >4 *2- Y < 4) ♦ 
142>*  <IZ~IX)4Y  <E>  4Y  <4>-Y(E>  4IR4  WR)/IY 

YF(6)  =  (RVR2*A2,T6*  (T1-T2)*YT*C  BT-(T14T2>  * XT4S3T 4 IXZ4  (YP(  4)  -Y<5>  4Y( 
16))4CIX-IY)*Y  (4)  4Y  (5>4Y<5>  4IR4WR)/IZ 
2MHT1-HT2) *XNT/IZ 

YF  ( .* )  =  Y  <1 )  4S  TFE-Y ( 2)  4CTHE4  SFE-Y  (3)  4CTHE4CFE 
Yf <8)=Y<5) 4CFE-Y(6)4$FE 
YP<9)=  <Y  (5  )  4  SFE4-Y  (  6)  4CFE)/CTHE 
YF(10)=Y  (4)4  (YF(9) 4STHE) 

YP(11)=Y  (1)4CTHE4CPSI4Y{2)  4  (SFE  4ST  HE  4CPSI- C  FE4  SFS  I )  ♦  Y  (3)  4  <  C  FE  4S  1 HE 
14CPSI4SFE4SFSI) 

YP(12)=Y(1)4CTHE4SPSI4Y  (2)  4  (SFE 4ST  F€  4SPS 1+  CFE4  CFS I )  ♦  Y  (3)  4(  CFE4STHE 
14SPSI-SFE4CFSI> 

444  CALCUL ATI CN  OF  EFFECTIVE  CM  AMO  CL 

GE  =  CL+(  (T1-T2)4YT4SAT  -  (T14  T2)4ZT4SBT  )/(RVR24A2) 

GC=CN+  ((Tl-T2)4YT4C8T-  (T 14-  T2) 4  X  T4S  BT  4  ( W  T  1-HT2 )  4  XH  T )  /  (RVR24  A  2 ) 

444  CALCUL  AT  3  CN  OF  FIRST  DERIVATIVES  CF  ALPHA,  BETA,  ANO  VR 

VCT=  (Y  (1)  4YP  (1)4Y(2)4YF(2)  4Y(3)  4YP(3)>/VR 
ALOT=  (Y(1)4YF(2)  -Y  (3)4YF(1))/{Y  (1)  44  2  +  Y  C  3 )  4  4  2 ) 

BCOT=  ( YP (  2  )  4  V  F  ~Y ( ? ) 4  V  C  T ) /O  L MX 
OLMX  =  VR4SCFT  ( VR4  4 2 -Y  (  2  )  4 42  ) 

RETURN 

END 
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SUBROUTINE  CLTPLOT 


THIS  SU  CRCUT I KE  FLCTS  THE  DESIRED  OUTPUT  AS  A 
FUNCTION  OF  TIME. 


*»#»##*¥*###** 


COMMCN/CMIN/EECMO,AV,S(E)  ,TRI HCPT ,S K ( 2 6 )  ,  ST  ART  ,  A  A ,  Ee ,  VOPT 
COMMON// AL  F  H  A  (40  2>,PETA(  40  2),VEL(4  02  ),P(402),LCCF,a(40  2), 

1R  ( 40  2 )  ,  ALT  (4  G2),  THETA  (  40  2)  ,PHI  (  402  ),  FSI  (  402)  ,Y0(7>  , TIME (40 2)  , 

2D  1ST (4C2) ?  TRf V (402) , AX  (4C2 ) ,Af (402), AZ (402) ,CEIE  (402  ), CEL A  (4C2)  , 
3DELR  (4C2i  ,  TGV  (40?)  ,THRST  (4  02) ,  CH"(  40  2)  ,CLE(402>  ,CNBEFF(4  02  )  , 
4LCF,CALF(4C2> , .  LOOK  4  0  2)  ,aEOOT  (  402) 

VV=VOFT*10C  .  *TRIMOFTMCCOQ  CO.*  BECMO 

CALL  PLOT ( C .  ,-12., -3) 

CALL  FACTOR  (  .c> 

CALL  SCALE (T 1NE* 10., LOOP, 1) 

CALL  SCALE  < A LFPA , 2 ., LCCF  ,1  ) 

CALL  SCALE <EETA, 2.  ,LCCF,1) 

CALL  SCALE (VEL,2.  , LOOP  ,1) 

CALL  SCALE (ALCCT, 2 ., LOOP, 1) 

CALL  SCALE  (EECCT, 2., LOOP, 1) 

CALL  SCALE  (P  ,2., LOOP  ,  1 ) 

CALL  SCALE  (C  ,2.,LOCP,l) 

CALL  SCALE  (R  ,2., LOOP,  1  ) 

CALL  SCALE(T»-ETA  ,2.,LCCP,1) 

CAuL  SCALE  (PH, 2.,  LOOP, 1) 

CALL  SCALE  (CELE,  2.  »LCCF  ,  1) 

CALL  SCALE  (CElR,  2.  ,LCCFf  1) 

CALL  SCALE  (CELA,  2.  ,LCCF,  1) 

CALL  SCALE  (A>,2.  ,LCOP,l) 

CALL  SCALE  (AK,2.  ,LCOP,i) 

CALL  SCALE  (A2,2.  ,LCOP,l) 

CALL  SCALE (C 1ST, 10  »,LCCF  ,1) 

CALL  SCALE (ALT  ,2. ,LCCC,1) 

CALL  SCALE  (TRAV,  2.  , LOOP,  1) 

CALL  SCALE (T FRST , 2 ., LCCF  ,1 ) 

L=LOCF+  1 
M=LCCP*2 

PHI  <L)=TC-V  (L  )  =  -190  . 

PH<M)  =  TGV  (H  )  =  13 0 . 

CALL  FLCT  (  (.  ,22. ,2) 

CALL  FLCT  (17  .,22. , 2) 

CALL  "LOT  ( 1 7  . ,  0.  ,  2  ) 

CALL  FuCT  (  C,  ,C., 2) 

CALL  PLOT (3.  ,3., -7) 

CALL  FLCT  (  C  .  ,  17. , 2 ) 

CALL  I  L  C  T  <  12. ,17.,  ?) 

CALL  FLCT  (12., C.  ,  2) 
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CALI  FLCT ( C •  ,0., 2) 

CALL  SYHBOL(?.,.25,.15,10HFUN  NUMBER  »  D .,  10 ) 

CALL  NUMBER  (9.5,  .2  5,  .15, VV  ,0.;  -1) 

CALL  FLCTC'.  ,12. 5, -3) 

CALL  AXIS <t,  ,C., lOHTI^S  -  SEC, - 10 , 10 . , 0 . , TIME ( L ) , TI f E CM) ) 
CALL  AXIS(C.,C., 11  HALFHA  -  DEG , 11 , 2.  , 90 .  ,  ALPHA ( L) , A LPH A ( M)  ) 
CALL  LINE  HUE,  ALP  HA,  LCCF,  1,0,  7  5) 

CALL  FLCT  (  C  •  ,-2.,-3) 

CALL  AXIS  ( C  •  ,  C .,  iOHTINE  -  SFC, -  10, 10 . , 0  .  ,T  IME ( L )  ,TI HE C M) ) 
CALL  AXIS  (C.  ,C., 19HALFHA-0CT  -  CEG/SEC, 19, 2. ,9 0.  ,ALCOT (L) , 
1ALD0T  (*-!)) 

CALL  LINE  (T  I  HE, ALOCT,LCCF,  1,0,7  5) 

CALL  FL CT  (  0.  ,-3.  ,-  3) 

CALL  AXIS  (  C  .  ,C.,10HTIME  -  SEC,  -  10,  10  . ,  0  .  ,T1ME  (L)  ,TI  ME  CM)  ) 
CALL  AXISCO.  ,C.,  10HBETA  -  CEG,  1  0,2  . »  9  0  .  ,  G  FT  A  (L  )  ,  BET  A  CM )  > 
IFCSTART.EC.  1.)  GO  TC  1 

CALL  SYM3CL  (  t  .2, 2. C, .15,12HPETA  CO  MM  AND  ,  0 , 0 , 12 ) 

CALL  NUMBER  <  C. 2,1. 7,  .  15,P£CMQ,  0  .,-1) 

CALL  SYM.BCLC  C.6,1. 7,  .  1 5 , 7H CEGR E ES,  0.  0,7) 

GO  TO  3 
1  CONTINUE 

CALL  SYMECIC  .2,1.7,.15,19HEETA  COMMAND  VARIES, C. , 19) 

3  CONTINUE 

CALL  LINE  CTIHE, BETA, LCCF,1  ,0,75) 

CALL  FL  CT  <  C.  ,-3.  ,-3) 

CALL  AXIS  <  C . ,C.,  lOHTIMf  -  SEC , -  10, 10 . , 0 .  , T IME ( L ) , TI HE (M) ) 
CALL  AXISCC.  ,C«,19H  BETA-DCT  -  CEG/S  EC  ,  1  9 , 2 . 0 , 9  0  . ,  0  E  COT  CL)  , 
1BECOT  (M) ) 

CALL  LINECTIHE, TED CT.LCCF, 1,0,7  5) 

CALL  FLCTCC.,-3.,-3) 

CALL  AXIS  (C.  ,C.,  10  HTIME  -  SFC,  -  10,  10  .  ,  0  .  ,T  IME  C  L  )  ,  TI  HE  ( M)  ) 
CALL  AXISCC.  ,C.,11HP  -  CEG/SEC  ,  11,  2.  ,  90  .  ,F  (L)  ,  P  CM) ) 

CALL  LINEUIHE,P,LCOP,1,0,75) 

CALL  FLCTC15  .,-4.5  ,-3) 

CALL  FLCTC  C.  ,22. ,2  ) 

CALL  FLCTC17  .,22.,  2) 

CALL  FLCT  (17  .,C. ,2) 

CALL  FLCTC  C.  ,C., 2) 

CALL  FLCT ( 2  .  ,2., -3) 

CALL  FLCTCC.  ,17.  ,2) 

CALI  FLCT  ( 1  2  . ,  17  , ,  2) 

CALL  FLCT (12  *,C. ,2) 

CALL  FLCTCC.  ,C.,  2) 

CALL  SYMBOL  (  7  • ,. 25,. 15 , 1 CH  RUN  NUMBER , 0. ,  10) 

CALL  NUMBERCE.5,  .2  5,  .15,VV,0.,*1) 

CALL  FLCTC 1.  ,13. 5, -3) 

CALL  AXIS (C,  ,C., 10HTIHE  -  SFC, -  10, 10 . , 0 .  ,T IME (L ) ,TINE CM) ) 
CALL  AXISC  C.  ,  C  • ,  11HO  -  DEG/SEC  ,  11,2.  ,90.  ,  CCD,  C(M)> 

CALL  LIKE  ( T  I  H E  ,Q ,  L  COF  ,  1,  C,  75) 

CALL  FLCTCC.  ,-3.,-3) 

CALL  AXISCC.  ,0*, 10  HT I *£  ■  SEC, - 10 , 10 . , 0 . , T IME CL ) , TI ME ( M) ) 
CALL  AXISC C.  ,C., 11  HR  -  CEG /SEC  ,  11,  2.  ,90  .  ,R  (U  ,  R  (M) ) 

CALL  LINE  (TIME, 7,1  COP  ,1»G»  75) 

CALL  FLCTCC.  ,-3. ,  -  3) 

CALL  AXIS CC.  ,C.,10FTINC  -  SEC,  -  10, 10 . , 0 . » T IME (L ) ,TI ME  CM) ) 
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CALL  AXIS  <  C  .  ,C.?  11HTHETA  -  DEG  ,  11,  2.  ,  90  .  ,  THET  A  (L )  ♦  THETA  ( M)  ) 
CALL  LINE  (TIME, THETA  ,  L  CCF,  1,0,75) 

CALL  FLCT(  0.  ,-2.  ,-3) 

CALL  AXIS ( C .  ,C., 10HTINE  -  EEC, - 10, 10 . , 0 .  , TIM? <L )  , TI NE ( M> ) 
CALL  AXIS  (  C  •  ,  0  • ,  9  H  PH  I  -  CE  G  ,9,  2  . ,  9  0.  ,  PH  I  (L )  ,  PHI  ( M)  ) 

CALL  FLCT  <  C •  ,  1.0, -3) 

0027J  =  1 ,LC  CP 
NC=2 


IF (J. EC. 1) GO  TO  26 
IF(APS(PHI  {J-1)-PHI{  J)  )  . 
XX=TIME(J) /T 1ME<M) 

TY=PHI  ( J )  / F HI  (V) 

CALL  FLCTOX  ,YY,NO) 

CALL  FLCT  (  C  .  ,-l.  0,-3) 
CALL  FLCT  (  C  .  ,-2.  ,  -  3) 
rail  ivTcif.  .r..int-TTMr 


GT  .50. )  N0=  3 


CALL  FLCT  (  C  .  ,-2.  ,- 3) 

CALL  AXIS  C  C .  ,C.,  10  FT  THE  -  SEC,  -  1 C,  10  . ,  0  .  ,  T  IME  ( L  )  ,TINE  ( M) ) 
CALL  AXIS  (0.  ,0.,  19HTUCNS  A  FT  GRAV  VE  CT  ,  19 , 2  . , 90  .  , TG  V  (L > ,  T 
CALL  FLCT  (  0.  ,1.0, -3) 

DC  13  J=1,LCCF 


CALL  FLCT  (  0.  ,1.0, -3) 

DC  13  J=1,LCCF 
N0=2 

IF(J.EC.l)  GC  TO  12 
IF(AB$(TGV  (J-l)-TGV(J)  )  ,GE. 

12  XX=TIME(  J)  /T  ])?(  M) 

YY=TG V ( J) /TG  V {►) 

13  CALL  FLCT  (  XX  »YY,  MO) 

CALL  FLCT((. ,-1.0,  -3) 

CALL  FLCT  ( 1 5  , , -4  •  5  ,-  3  ) 

CALL  FLCT ( 0 .  .22. .2) 


TG  V ( M  ) ) 


50.)  NO=3 


CALL  LLCT ( lb., “4 .5 , 
CALL  FLCT ( 0 .  ,22. ,2) 
CALL  FL CT ( 17 , , 2? • , 2) 
CALL  FLCT(  17  .,0. , 2) 
CALL  FLCT (C.  ,0., 2) 
CALL  FLCT  (  2 .  ,2., -3) 
CALL  PLOT  <  C •  ,17. ,2) 
CALL  FLCT  ( 1 2  • ,  17 ,  ,  2) 
CALL  FLCT  ( 12.,  0.  ,2) 
CALL  FLC  T  {  I ,  ,  C.,  2) 

At..  a  i.  >  i  n  a  i  #  •»  a  p 


CALL  *•  LL  T  (  .  ,  U,  ,Z  » 

CALL  PL C T  C  C .  ,C., 2) 

CALL  SYM9CL(7.,.  25  ,.  15,  1  CH  RUN  NUMBER  ,  0 .,  10  ) 
CALL  NUMBER  («.5,  .2  5, ,  15,  W  ,0.,  -  1) 

CALL  PLCT  ( 1 ,  .12.5,  -3) 

CALL  CXTS(f  .  19HTTMP  -  CRC .  -  1 1 .  1 0  .  .  0  .  .  T  I 


V  rui.  •  w  I  *  *  ■  .12.5,  — 3) 

CALL  AXIS(C.  ,t., 10HTIMC  -  SEC, -  10, 13 . ,0 . , TIME (L )  , TI ME (M) ) 
CALL  AXIS(G.  ,C., 1 9  PEL  E  \,  A  TO  K  9ISF  -  DEG,iS,2.,90«,OELEiL) ,9 
CALL  LINE  CTIf'E,nELF,LCCP,l  ,0,7  5) 

CALL  FLCT  (  C.  ,-2.  ,  -7i 

CmlL  AXTM  C.  ,C.,  10HT  IuE  -  SEC,  -  10 ,  10  .  ,  0  .  ,  T  I  ME  ( L  )  ,  T I  PE  i  M)  ) 

CALL  AXIS < C.  ,C.,  I7r'!  CCER  flSP  -  OEG , 17  ,  Z  .  ,  9 C . , C £ LR  (L ) , DEL 

CALL  LINE  (1IRE,3ELF,LCCr.i  ,0,75) 

CALL  FLCT  (  C.  ,-2.,-3) 

CALL  AX  IS  {  f  .  ,C.,  1 0  t-T  I  v  E  -  SEC,  -  1 0,  J  0  .  ,  0  .  ,  T  IMF  ( L  )  ,  TI  ME  ( M)  ) 

CALL  AXIS(C.»C,,  1 3  HA  I L  F  R  CM  OISP  -  CEG  ,1  8  ,2.  ,90  .  ,  CEL  A  (L  I ,  DE 

CALL  LINE(TIME,')ELA,LCCF,1  ,0,7  5) 

CALL  PICT (0.  ,-2. 7) 

CALL  A  X IS ( G «  ,C.,  10(TTME  -  SLC,  -  10 ,  10  .  ,  0  .  ,  T I  ME  ( L )  ,  T I  ME  ( *>  ) 
CALL  AXIS(C.  ,C.,  17K  FlCriTV  -  FT/SEC  ,  17 , 2.  ,90.  ,  VEUU  ,  VEL 
CALL  LINE  (T  IMF. ,'/ EL  ,L  CC  F  ,  1,  0,75) 

CALL  FLCT  (  C  .  2.  .  -2) 


LR(I')) 


L  A  (  M  5  ) 
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CALL  AX  I S ( 0.  ,C., 13H0ISTANCE  -  FT,- 13 ,10 . , C. ,OIST (L) ,OIST (M ) ) 
CALL  AXIS  (C.  ,C.,  13HALTITUDE  -  FT,1 3,  2. ,  9  0.  ,  ALT  ( L  >  ,  ALT  <M)  ) 

CALL  LINE(CIST, ALT, LCCF,1, 0,75) 

CALL  FLCT(15.,-4.5,-3) 

CALL  FLCT<  0.  ,22. ,2) 

CALL  FLCT  ( 17  •  ,22  .  ,  2) 

CALL  FLCT ( 17  •  ,0. , 2 ) 

CALL  FLCT ( 0  «  ,C.,2) 

CALL  FLCT (3.  ,3., -3 ) 

CALL  FLCT  <  C-  •  ,17.  .2) 

CALL  PLCT(12  .,17.,  2) 

CALL  FLCT (12. ,0., 2) 

CALL  FL CT  (  0  .  ,C., 2) 

CALL  SYM9CL(7.,.25,.15,1CHRUN  NUMBER  ,  0.  ,  10) 

CALL  NUMEER  («.5,  .2  5, .  15,  VV  ,0.»  -  1) 

CALL  FLCT(1«  ,13. 5, -3) 

CALL  AX  IS ( C •  ,  C ., i 0  HT IM  F  -  SEC,  -  10,  10 . , 0 . , T IMF ( L ) ,TINE ( NS ) 

CALL  AXIS  (0.  ,C*,  13HX  ACCELEPAT  ICN  -  G,  1  8 , 2  .  ,  SO  .  ,  A X  (L)  ,  AX  (M  ) ) 
CALL  LINE (TIME, AX, LCCF  ,  1,0  ,75) 

CALL  FLCT  (C.  ,-2.  ,-3) 

CALL  AXIS  <  C  •  ,C.,  13HTIME  -  SEC,  -  10 ,  10  . ,  0  . ,  TIME  ( L )  , TI  NE  ( »)  ) 

CALL  AXIS  (C.  ,C.,  16MT  ACCELERATION  -  G, 16, 2., SO., AY  (l),  AY  (M)) 
CALL  LINE(TUE,AY,LOCC,1,0  ,75) 

CALL  FLCT  (  0  .  ,-3.  ,  -  3) 

CALL  AXIS (C.  ,C., 1 0  HT I v  c  -  SEC, - 10, 10 . , 0 . , T INE ( L )  ,T I NE  ( N) > 

CALL  AXIS  <  C.  ,C.,  16HZ  ACCELERATION  -  G,  1  fl  ,2.  ,90  .  ,AZ  (L)  ,  AZ  (1  )  ) 
CALL  LTNE(TINE,AZ,LOCF,1,0  ,75) 

CALL  FLCT  (C.  ,-3. , -2) 

CALL  A  XIS  (  C  •  ,C.,  10HTIVF  -  SEC,  -  10,  10  . ,  0  .  ,  T  I  ME  (L )  ,  TI  ME  ( M)  ) 

CALL  AX  IS  C  C  .  ,C.,  12NTNRLST  -  L9  S  ,12  ,2  . ,  9  0  .  ,  THRST  ( L)  ,  THRST  (M  )  ) 
CALL  LINE(TINE,THRST,LCOP,  1,0,7  5) 

CALL  FLCT  (  C  .  ,-3.,-3) 

CALL  AXIS(C.,C»,i3H0ISTANCE‘  -  FT ,- 13  ,1 0  .  ,  C  . ,  OIST  (L)  ,OIST  (M  ) ) 
CALL  AXIS  (C.  ,C.,  l&HCRCSS  RANGE  -  FT,  16 , 2.  ,  90.  ,  TRAV  <L)  ,TRA  V  (M  ) ) 
CALL  LINE  (DIET, TRAV,LCCP,1  *0,75) 

IF(START.NE.  1.)  GO  TC  5 
CALL  SCALE  (C  ALF,  9.  ,LOF,l) 

CALL  SCALc  CC  N  E  ,4.  ,LCF,  1) 

CALL  SCALF  (CIE  , 4.  ,LOF,  1) 

CALL  SCALE(ChEErF,4.  ,LCF  ,1) 

CALL  FLCT  (15  .,-4.5 , -3) 

CALL  PLCT ( 0 •  ,22.  ,2) 

CALL  FLCT  ( 17  , ,  22  •  ,  2) 

CALL  FLCT  (17. ,Q.  ,2) 

CALL  FL CT  (  G •  ,G., 2) 

CALL  FLCT  (2.  ,3., -3) 

CALL  FLCT ( 0 •  ,17. ,2) 

CALL  FLCT  (12. ,17,,  2) 

CALL  FL  C  T  ( 1 2  , ,  C.  ,2) 

CALL  FLCT  (C.  , C.,  7) 

CALL  S  Yf-GCL  <  7 .  ,.25,.  15, 10HRUN  N'L'Mn  ER  ,  0.  ,  1 C  ) 

CALL  NLMBERP-.5,  .25,  .  15 , VV  ,0.,-  1) 

CALL  FL  C  T  (  2  .  ,12., -3) 

L=LOP+l 
H=L  + 1 
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CALL  AX  IS ( 0 « « C . , 1 1 HALFH A  -  PEG  , -11 ,9  . ,  0  . ,  CALF  (L)  ,  CALF  (  K>  ) 

CALL  AXIS(  0.  ,0.,  HHCN-oETA -EFF  ,  11,  4.  ,90  .  »  CNE (L ) »CNE (M>) 

CALL  LINC(CALF,ONE  ,LCP,  1,-1, 3) 

CALL  FLCT  (  C  •  ,-5.,-3> 

CALL  AXIS (0.  ,C., tlHALFFA  -  DEG  , -11 ,9  . ,  0  .  ,  CALF  (L),  CALF  m> 

CALL  AXIS ( 0.  iun  11FCL-PETA-EFF  ,11,4*  ,  90  •  ,  CLE  (L)  ,CLE  (M)) 

CALL  LINE(CALF,CLE,LCF,1,-1,3> 

CALL  FLCT ( 0.  ,-5.  ,-3> 

CALL  AXIS ( C •  » C«»  11  HALF  HA  -  DEG  , -11 ,9  . ,  0  .  ,  CALF  (L),  CALF  (F>> 

CALL  A  XI S  C  C  •  ,  C.,  15FCN-PE TA -DYN- EFF ,1 5 , 4 .  ,90.»CNEEFF(L),CN9EFF(P>> 
CALL  LINE!  CALF  ,CNBEFF,LOF  ,1,-1, 3) 

5  CCNTIMJE 

CALL  FACTOR  (  1.) 

CALL  FLCT (12  .,-12.  , -3) 

RETURN 

END 
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SUBROUTINE  S  FEVAL ( Y, YP ) 


***  THIS  SU8  ROUT IN!  PEFFCRMS  THREE  FUNCTIONS  t 

1)  C  AR  IS  C  h  CF  THE  SIGNS  OF  THE  FIRST  ANC  SECCMC 

CER1VATIVGS  OF  BETA  FOR  OET  E  RH INAT I C  N  CF  THE 
CEFATTORE  FCINT 

2)  CHECK'S  THE  VALLE  OF  PSI  AS  A  SECCNC  CEPARTURE 
CRITERIA 

3)  CALCULATION  CF  AN  EFFECTIVE  C-N’-QET  A- DYNAMIC 


OIMENSION  Y  ( 12),  Y®  H2> 

CCMMCN/CNAIN  /  FECMD , A  V , A VC5  , AVO T , JA ,A C, PI  ,TCPI , TRIHCr  < , NE AR  , 
1THPUS1,TH0LS  c  ,A1,A  2,  A  3,  El,  e2,9  3,Tl  ,T  2 ,  F  EC  ,  S  IGHC ,  THE  C,  ANGLE  , 
2GA,G3,GC,OE,CA,OR, T°T 1 ,TRT  2 ,AL , PE, AL PHA T , ST  ART , CNEE ,CLSE ,V CPT , 
6SFE,CFE, VR ,CYCR» CNCRf CHCF, THRUS ,TP IM , BET  AT , WTl , HT2, XWT,CON  T ,ZHT  , 
7ALOT,PECT,ALCFT,3ECPT ,TCFT ,XT,  YT,ZT,  IX , I Y , I Z ,1 XZ , IR , WR ,M AS S , G  ,  S  , 
8B,C,CX,CXOE,C7,CZOE,CW.CMOE,CMQ,RVP2  ,ROH,STHE,CTHE,TVGAIN 
CCMMCN//ALFHA  (AC  2)  ,6ETA(<.O2),VEL:4  02  ),P(4C2)fLCCF,C(4O2), 
lF(402),ALT(4C2i,THEfA(4C2)  ,®HI  (  402 )  ,  PSI  (  4  C  2  )  ,  YO  ( 7  )  ,  TIME  (  40  2 )  , 
20IST  (4  02) ,TRtV'402>, AX  (402  ) ,AY (402) , AZ(  402) ,CELE (40  2) ,CELA  (40  2) , 
30ELR(4C2)  ,TG  W43  2)  ,THPST  (4C2) ,  CNE(  40  2)  ,  CLE  (4G2)  ,CN9EFF  (402  )  , 
4t.CP,CALF  (4G2  )  ,AL 001(4  02)  ,PEDOT  (402) 

REAL  IX,IY,I2,IXZ,IR,HASS 

no  FORMAT (1H1, 7 >,*J*,  12X  ,  *C EL S* ,1 IX,* OELCN*  ,  1 IX  ,* C ELCL* , 3X, *C N5E FF « , 
14X,*ALCHA*,E>,**)ETA*,//) 

EECM=EECMC  *  A  V 

10  0  FORMAT  (1H 1,7 >,*J*, 12 X , *CEL e*, 1 1 X,* OELCN *  ,  1 1 X ,*DELCL *  ,  3X, *C NBE FF*  , 
14X,*ALFHA* ,5>,*9ETA*  ,//) 

101  FORMAT  (4X,I4,2(4X,E12.5),3  (4X,  F5.2>> 

103  FORMAT  < 1 H 1 ,////, 4X ,*  THE  FO INT  OF  DEPARTURE  OCCURS  A7*,//,12X, 

1*  ALPHA  s  *,F10.5,*  DEG  REES  *  ,//  ,8  X ,  *  ANC  BETA  =  *,F10.5, 

2*  OEGREES*  ,//,9X,*AT  T IM E  =  *, 3X , FI  0 , 5  ,  *  SECONDS*, 

?//,5X  ,+ ANC  AT  P£TA  CCHMANT)  =  *,F10.5,*  CEGRGES*  ) 

104  FORMAT  (1H1  ,////,  *  NO  DEPARTURE  DEFI  NEC*) 

105  FCRHAT ( //// ,  *  BETA /° E7 A- COT  CONDITION*) 

106  FORMAT  (////,  *  PSI  CONDITION*) 

HAX=LCOF-l 

1=0 

00  2  J= 1 ,M AX 

IF  (TINE  ( J )  . L  T  .5. >  GO  TC  1 
BFOOT  =  EECCT  ( v  > -9EOCT  ( J-l ) 

TEST  =  GECCT  (J)*EEOOT 
IF (TEST  »GT  .  C  •)  GO  TC  20 

1  CONTINUE 

2  CONTINUE 
PRINT  104 

20  1=1+1 
JrJ+l 

0EDOT  =  EcOOT  (JH9F.r>0T  ( J  + 1 ) 

TEST  =  EECCT  ( J ) *?E  00  T 
IF (I. EC.  3)  GC  TO  3 
IF(J.cO.MAX)  GO  TO  1 
IF(TCST.GT.C.)  GO  TC  20 
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1=0 

GO  TO  1 

3  PPINT  103,AIFM(  J)  ,RETA(J>  ,TIME(J>  ,BECMC 
PRINT  105 

4  00  5  J=2 ,LCCF 

TEST= A8S  CP S I  <J)*36C. ) 

IFCTEST. GT . 1 5  .)  GO  TO  6 

5  CONTINUE 
PRINT  104 
GO  TO  <3 

6  PRINT  1C3,AIFFA(  J)  ,BETMJ)  ,TIHECJ>  ,BECMO 
PRINT  106 

9  CONTINUE 
PRINT  100 
MAX= J-5 

0FLT=TINE ( 3) -TIME ( 1) 

1=1 

00  8  J  =  1 , M A  X 

IF (GET A  (J)  .G1.5. )  GO  TC  7 
IF  CPFT  A ( J)  .EC.O. )  GO  TC  7 
CL8=CLE(-J)  /BETfi(  J) 

CNR=CNE(J)/EETA(J) 

AIF  =  ALFHQ  (J  )  /AV 

CNBEFF(J)=  ( C  C*C  CS  (AlF)  -  (I  7/IX  )  *CL  Br  SIN  CALF)  ) 

IFCAGS  (CNEEFF  (J)  ).  CT.  0.10)  GO  TC  10 
CKBEFF (I)=CNFEFF (J) 

CLE(Ii=CLn 
CN'E  (I )  =CNB 
CALF(1)=ALF*AV 
1=1*1 

7  CONTINUE 

8  CONTINUE 
L  CF= I- 1 
RETUPN 

10  PRINT  101,  J,  EECOT  (  J)  ,  ONE  (J  )  ,CL  E  CJ>  ,CN9EFF  (J)  ,  ALPHA  (J) ,  BETA  (J) 
GO  TC  7 
END 
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SUBROUTINE  STCEE(Y,YP) 


DATA  STORAGE  FOR  PLOTTING  SUBROUTINE 


OIHENSICN  Y ( 1 2 ) »  YP  (1 2 ) 

COMMON/ CM A  IN  /  EEC MO ,A  V  *  AVC5 , AVOT , JA , A  0, P I  ,T CPI , Tfi IMCFT  ,  NE AR  , 
lTHRUSl,THRUS2,ai, A  2, A?,ei, E?,3 3 ,T1 ,T  2,R EO , SIGHC , THE C , ANGLE , 
2GA,GB,GC,nE,CA,0R,TRTl  ,TRT2,Al  ,  Er  ,  AL  PHAT  , START  ,  CNEE  ,CLEE,VCPT, 
6SFE»CFE,VR,CYC‘J,CNrR,0HCP,  THRUS  ,TR  IM  ,  0£  T  AT  ,  WT1 ,  HT2,  XWT  ,C  CN  T  ,  Z  NT  , 
7AIDT,EECT,41CPT,  3E CPT , TOFT ,XT,  YT ,Z T,  IX , I Y , IZ ,IXZ , I R ,WR ,M AS S , G  , S  , 
8e,CtCX,CXCE,C2,C70E,C*',RN0E,CM0,RVR2,»0H,STHE,CTPE,TVGAIN 
COMMON/ /ALCH  A  (A3  2)  ,O£7A(4G2),VEL(4  0?),P(4O2),LCCF,Q(4O2), 

1R(4  02),ALT(4(2),THETA<4C2),phI(  402),PSI(4C2),YC(7),TIME(  40  2>, 
2DIST  (402)  ,TRAV(402),  AX  (4C2  )  ,A.Y  (402  ),  AZ(  4  02  )  ,CELE(40  2)  ,CELA  (4  0  2  )  , 
30FLR  (4  02)  ,TGV  (40  2)  ,THFST  (402)  ,  CNE(  40  2)  ,  OLE  (4  02  ) ,CN3EFF  (4  02  )  , 

4L CP, CALF  (4C2  )  ,  AL  00  T  ( *  0  2 )  ,9  EOOT  (402) 

REAL  IX,IY,I2,IXZ,  IR,NAS? 

J=GA 

AX(J)  =  (  (YP  (1  )-Y(2)  *Y  (6)  +  Y(  3)*Y  (5))  /G)  kSTHE 
AY(J) a((YP(2)-Y(3) *Y (4)*Y( 1)*Y (6))  /G)-CThE*SFE 
A7(J)*((VP(2)-V(i)  *Y(5)+Y(2)*Y  (4))  /G  )  -C  THE*CFE 
CELE  < J ) -CE 
OELR ( J) =CP 
DEL  A ( J) =CA 
\  HRST ( J) =  71+12 

IF(Y(C)  .GTrPI)  Y(9)  =  Y(B)*-TCPI 
IF(Y<c).LT.-FI)Y  (9)  =  Y(c>  +TCPI 
TGV(J)=Y(9) *AV 
OIST  ( J) =Y ( 1 1 ) 

TRAV ( J) =Y ( 1 2  ) 

ALPHA ( J) =AL* AV 
ALOOT  (J) =ALCT*AV 
BETA  (J)-PE*A  V, 

8  FOOT  (J)”9EC1#A\/ 

VEL(J)=VF 
P(JI  =  Y  (4)*A V 
0 ( J) =  Y (5)* A V 
R ( J )  =Y  (6)*AV 
ALT ( J) =Y  (7 ) 

THETA (J) -Y (5  )*A7 

IF(Y(10).GT.FI)  Y(  10)  =  Y  (10  l-TOPI 
IFCY(IO)  .L7. -FI)  Y(  IQ)  *Y(10  ) +TOPI 
PH(J)=Y(1C)  7 A V 
PSI(J)=Y(9)*AV/360. 

CLE ( J) =GB 
CN'E  ( J)  =GC 
RETURN 
END 
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SUBROUTINE  T  FIH1 <Y  ,YF  > 


***  AIRCPAFT  IS  TRINMEO  FOR  STRAIGHT  AND  LEVEL,  0 LIMPING, 
OR  OESCE  NCING  FLIGHT 


OIMENSICN  Y  ( 12), Y0  (12) 

CCMMCN/CMAIN /EEC  MO,  A V, AVC5 ,AVOT ,JA , A  0, P I ,TCPI , T R IMOFT , NE AR  , 
1THRUS1,TFRLSZ  ,A1,A2,  A2 ,91,  E2,9  3  , Tl  ,T  2  ,  F  EC  ,  S IC-HC  ,  THE  C  ,  ANGLE  , 
2GA,GP,GC,0E,  CA,DR,  TRT1,ThT2,AL  ,  nE,  aiPHAT,START,CNEE  ,CLEE,VCPT  , 
6SFE,CFE,  VR,C  YCR,CNDR,CFCF,  THRUS  ,T=>  IM  ,  PET  AT  ,  KT 1 ,  W T2 ,  X WT  ,CCN  T  ,  Z  M T  , 
7AIDT ,EECT ,ALCFT, PECPT ,TOFT ,XT,  YT,ZT, IX, IY, IZ,IXZ,IR ,WR,MASS,C-  ,S  , 
80,C,CX,CXCE,CZ,CZDE,CM,CNDE,CMQ,RVR?  ,R0 H  , STHE , C T ME , T VG AI N 
CCMMCN/CCCEF/FCY (2 1 , 5 ) , E CN (21 ,  9 > ,E CL (21 , 9 ) , ECM ( 2 1 , 9 ) , EOX (2 1 ) , 
1ECXCE  (21) ,£C  2CF(21 ) , 2 C  MC  E ( 21),  EC^n (2 1) , EC YCP (2 1 ) , ECN DR t  ' 1)  , EC LC F  (  2 
21)  ,ECYOA  (21)  ,  ECHO  A  (2  1 )  ,  ECL  CA  (2  1  >  ,E  CY  P  <  2  1 )  ,ECNP<  21)  ,EC'.F(  21  )  ,ECYR  (2 
21)  ,  ECN'F  (21),  ECLRC21)  ,FC?(2  1) 

RFAL  IX,  IY  ,  I  7  ,  IX  Z,  IP,FASS 

100  FORMAT  ( 1F1 ) 

101  FORMAT  ( 2 2H  IE  ITERATIONS  ON  TRIM,/ IX  ,7HUCCT  =  ,E12.5,8H  WDCT  -  , 
lE12.5,fiH  CCCT  =  , E 12  *  E  ) 

102  FORMAT  (*  AIR  CF  AFT  IS  TRIMMED  FOR  STRAIGHT  ANC  LEVEL  OR  CLIPPING  FL 
1IGHT*,/,*  •,62(lH*>,////> 

103  FORMAT  (23M  TRIM  CZ  NOT  03  TAIN  ABLE ,/ 14H  DESIRED  CZ  =  ,E12.  5) 

10 A  FORMAT  (SI-  CE  »  ,FE,1> 

105  FORMAT  (  9F  CESCZ  =  ,  F7.3) 

106  FORMAT  (<?H  TFFLS  =  ,E12,5> 

108  FORMAT  (5F  ECZ(,I2,4F>  =  ,E12.5) 

109  FORMAT  (9H  YF(1)  =  ,  E 1 2 . 5,  ?H  Y  P  (3)  =  ,E12.E,5H  YP(5)  =  ,E12.5) 
PRINT  100 

PRINT  102 
T=0. 

AL=ATAN2  (Y  (3  ),Y(  1)  ) 

CCUNTEF=0. 

Y(4)=Y(5'  =  Y(6)-Y (9)  =  Y  (10)= Y(2) =0. 

7  CONTINUE 

COUNT  ER  =  CCUME  Ft  1, 

Y ( 8) =  A  L 

CALL  GYRATES  IT  ,Y,Y  F) 

CF=- ( CMt  (Tl  +  12>*ZT/(RVR2*A3*MASS>>  /CHDE 
HPITE  (6,104)  CE 
CALL  GYRATES  (T,Y,YF) 

OESCZ=-  (  (G’CCS(AL)  /<9VR2*A  1))  +  CZDE*DE) 

WRITE  (6,105)  CESC? 

2  CONTINUE 

IF(JA.LT.1.CR.JA.GT. 21)  GO  TO  5 
IF<ECZ(JAi .EC.CESC7)  CC  TO  3 

IF  (ECZ(JA).  CT.OESCZ.  AND.ECZ(JAtl)  .LT.OESCZ)  GC  TO  1 
JACC=C 

IF  (ECZ  (JA  )  •  CT  .3  ESC?)  JACO-1 
IF(ECZ(JA)  .LT.CESC7)  JA00=-1 
JA= JA* JACC 

WRITE  (6, 1 C8  )  JA , ECZ ( J A) 

GC  TO  2 
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CONTINUE 

fiLs <JA«3)/AVC5 


GO  TO  4 

a“T(  U»-3>  -  «C2.  JA  >-75507,  /.EC  7  «JA.l  l-ECZ  ( JA„  )  /MOS 


GO  TO  4 
CONTINUE 
WRITE  ( 6  *  1 C  3  ) 


CESC7 


CHOP=l. 

RETURN 

CONTINUE 

CVR^T  ( 1)  **  £  ♦  T  2 

y  (l)=SORT(CVR/(i.>TAN  (A 

Y  (3)=SCRT (CVR-Y(  1)  **  2  ) 

Y  ( 8  )  “  A  L 

THRUSMC-*$  IN  ( Al)  -RVR2 *  A  1¥C  *fCX  0  E*OE)  *MASS 
T1=TRI  1*THcU; 

T2=TRT2’THRLS 

WRITE  (6*106)  THRUS 
CALL  GYRATES  ( T  tY  » Y  c> 

.  AS  $  (YP  (3>)  .GT. 2 • »CR 

1  GO  TC  7 
RFTURN 

WRITE  (6*101 )  Y9(D»YR(3)»YP(5) 

RETURN 

PKn 


,ARSCYP(5)  ).GT  . 


( 
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SUPRCUTISE  7P1*2(Y,YP) 

OIKENSICN  Y ( 1  ? ) i Y°  ( 1 2 ) 

C  C  K  K  C  N  /  C  *•  A  I  f ,  /  i  E  S  M  0  ,  A  V  ,  A  V  C  5  ,  A  V  0  T  ,  J  4  ,  ft  0 ,  P  I  ,  T  C  F  I ,  T  R  I  M  C  F  T  ,  N  E  A  R  , 
1THRUS1,TH°US2,A1 , A  2, fl 1 , 5 1 ,  E2, 3 3 ,Tl ,T 2 , F E C , $ IG« C , THE C , ANGLE , 
2GA,GP,GC,DE,  [A  ,QR,  TRT  1  ,  TPT  2 ,  AL  »  Cc.,  M  i-HAT  , ST  ART,  CNEE  ,CltE,V  CRT  , 
6SFE,CFE,  VR,CYCR,CNCt>,OCP,  THRU  $  t  T  R  XM  ,  *E  T  AT  ,Ml,W2,m,CCNT,ZMT, 
7AL0T,PEDT,^LCFT,'5E  Ct>T,  TOFT  ,XT,  YT,ZT,  XX , XY, X7 ,XX2 , XR *WR,M AS S*G »S* 
0e»C»CX,CXOE  ,C2  ,<:7DF,CJ'  ,C^OE,CMQ,RVR?  iRCH  ,STHE,  CTPE,  TVGAIN 
CCMMCN/CCCEF/ECY (21, E)  ,ECN  (?1,  9),ECL  (21  ,  E)  ,  ECM  (  21 , 9  )  *  ECX (2  1)  , 
1ECX0F  (21)  ,EC2CE(?l),ECVOE(21),  CCMQ  (2  1)  ,  ECYCR  (2 1 )  ,  ECNOP  (2 1)  ,  EC  LC R  ( 2 
2l),ECY0A(21),ECN0A(21),ECLCA(21),tCYF(21),ECKP(21),ECLF(21),ECYF(2 
31>,ECNR(21),tCLR(21)  »CC7  (21) 

REAL  IX,IY,I2,IX2, 1° , M ASS 

♦  ♦♦♦♦♦ENTER  TUFN  TRIM  FERE  *♦*♦♦♦* 

RETURN 

ENO 


ooono 


■Wfn-I.  '■im  uiMumi.iiijpii 


GAM/AE/73-14 


SUORCUTIKC  FRINTl(Y) 


♦**  AIRCRAFT  F  ARA  METE  R  c  ANO  INITIAL,  TRIMMED  CONDITIONS  ARE  PRIMED 


DIMENSION  Y (  1 2» 

CCPPC^/CPAIN/EECM:i  ,AV,  AVC5  ,AtfOT  ,JA,A0,PI,7CPI,TRIM0FT,NEAR  , 

1THRUS 1,  T  HR  l  S  ?,  Al  ,  A  2,  A  3,  PI, E2,n 3 ,Tl ,T 2, FEC, SIGHC , THEC , ANGLE , 
26A,Ge,GC,CE,  CA,[»?.  TRT1  ,TRT2,AL  »PE,  Al.  °HA  T, START, CNEE,  CL  EE  ,VCP7  , 
6SFE,CFE,VR,C'iCR,CNrt?,CFCF,  THRU  S  ,TR  IM  ,  BE  T  AT  ,  WT 1 ,  RT2  ,  XHT  ,C  CN  T  ,  Z  N  T  , 
7ALOT,PEDT , ALCFT, S£ CPT , TCFT  ,XT#  Y T,Z T, IX , I Y, IZ , I X 2 , 1 R ,KR ,M AS S , G  , S  , 
8B,C,rx,CX0E ,  C7,CZ0F,O,CNDS,CMQ ,RVR2 ,  RO H  ,  S T HE , CTHE , TVGAIN 
REAL  IX, IY, 12, 1X2,  IR,NASS 
ALPE=AL*AV 

WRITE ( 6 , 1  Q  7 )  Y(l)  , Y ( 3 )  ,  ALP  E 
IFIMsl. 

PRINT  110, IX  ,IY,IZ ,IXZ,MASS,S,C ,9, IR ,WR ,T1  ,T2, THRUS 
PRIM  111,  (Y  (J),  J=l,6)  ,THEC,FEC,SIGHC,Y  (7)  ,CE,CA,CR 

11 0  FORMAT (1  HI ,////, 4< IX , 1 U0  (1 HM , /> ,4 (1 X,4H**** ,9  2X ,4H *♦♦*,/)  ,1X  ,4h** 

1**,30X,32H***  AIRCRAFT  PARAMETERS  •♦♦ , 3 CX , 4H** ♦♦ , /, 3( 1 X , <■ F ♦♦ * 
2*,92X,4H»»**  ,/),  IX  ,4H*»»»,  5X,*  IXs*  ,E14.  3 , 4X  ,  » I Y  s*  ,  E  14 . 8,  4X  *  IZ  =  *  ,  E 1 
34.8,4X,*IX7s*,El4.  8,6X,4l-****,/,iX,4H****,9?X,4H****l/,iXl4H****f5 
4X,*MASS**,E1  1.5, 4X,* WING  A REA= * ,E1 1. 5 , 4 X  ,»CWCRD =  ♦ , E 1 1 . 5, 4X , ♦ S F A  Ns* 
5,Eil.5,5X,4h****  IX,  4h****, 9  2  X,  4  H*  ***,/,  1X,4H****  ,5X,  NOTARY  1ER 

6MS*  INERT  I  As* , E 12.6, 4X,  *FR£  OUENCY =  * , E 12 .6 ,2 4X , 4H **♦* , /,  IX, 4F*** 
7*,92X,  4H****  ,/,lX,  Ml****  ,«?X,*SMGIME  THRUSTS  l  LEFT  EMINE  =  *,E12. 
86, 4X,  ♦RIGHT  E  NGI  \*E=*  ,  E  1 2 .6 , 15X  ,  4H*  ♦♦♦,/,  1  X  ,  4H*  *♦♦,  2  4X ,  *T  CT  AL  THRUS 
9T=*,E13.  7, 4  2>,4H****,  3  </,lX,4H****  ,9  2X  ,  4HJ  ♦♦♦)  ,3(/,lX,10  0(  1 H  ♦  ) )  ) 

111  rCRMAT  ( 

14  (IX,  4H****,S2X,4H**  ♦♦,/),  1X,4H****,  3  OX,  3  IF***  INITIAL  CCNCI1IC 
2NS  ***,3 IX  ,4F****,  /  ,3(1X  ,  4H*  **♦,  92 X, 4 H ****,/)  ,1X,  4H****,  EX  ,  *1  (C) 
3=*,E14.8,4>,*V<0)  =♦,£  14.8,  4X ,  ♦  W  ( 0)  **  ,E  14.8, ?2X , 4H** **  , /,  IX  ,  4H****, 
492X,4F****  ,/  ,1X,  4H****  ,5X,  *P  CO  )  =  ♦,  FI  4. 8 , 4X,*C  (C  ) -*,  El  4. 8 ,4  X,*R  (C)  = 
5*,£14. 8,22X, tH**** ,/ , IX, 4H ♦♦♦♦ , 92X ,4H** **,/,lX ,  4F****,5X ,♦ THETA  (0) 
6«*,cl4.  fi  ,4  > ,  ♦FHK0)s*,El4.  8,4X  ,  *PSI(  0)a*,El4.3,  14X,4H***V  ,1X  ,4H** 
7**,  92  X,  4F*  ***,/,  IX  ,4F****,  EX,*  INITIAL  A L TI TUCE  =  v , E 1 4 . 9 ,5 6X , 4H * ♦ **  , 
8/,lX,4H****,C2X, 4H ♦♦♦*,/, IX, 4H**** ,5 X,* ELEVATOR  QISFLACEMZ NTs*  ,  E14 
9. 6, 4X, ♦AiLER  CN  9 I*FL ACEMFN T  =  *, E 14.  S, 12X , 4H* ***,/, IX , 4H**** ,5 X  ;  ♦ RUG 
10EP  DISPLACE HFNT=* ,E14.8,5 3X,4H****,  /, 3 ( IX , 4H** ♦♦ , 9 2 X, 4H ** ♦♦,/) ,4< 
21X, 100(1 H*),  / ) : 1H 1 ) 

107  FORMAT  ( 5H  U  =  ,E12.5/,EH  W  =  ,E12.5/,5H  A  s  ,E12.5) 

RETURN 

END 
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1 

30 

2 

3 

4 


9 

10 


20 

21 

22 

23 

24 
11 


•AES  (H))  1,3,3 


K 

,F0) 


12 

13 

14 

15 
17 


16 

19 


xc=x 
x=x+ox 

W=Q.5*  (X-XG) 

H  =  H4-H 

IF  CAES  (X-XC) 

OC  4  1=1, N 
Y0(I)=Y(I) 

HT  =  H 
XT=X0 

00  5  1=1, N 
Y7 ( I ) = Y  0  (I) 

ASSIGN  e  7 C  K 
GO  TO  20 
CO  7  1=1, N 
YF(I)=Y(I> 

HT=0.5*H 
ASSIGN  9  TC  I 
GO  TO  20 
DC  10  1=1, N 
YT(I)  =  Y  (I) 

XT  =  X  0  +  F  T 
ASSIGN  11  TC 
CALL  F(XT,YT 
DO  21  1=1, N 

YU>  =  YT<IMC.5’MT*F0<I> 

CALL  F  (XT  +  C. E #NT , Y :P1 ) 

5ai^!u!*H'l.  207ia7ei*FO<I>  +.292  a93219«Fl<»> 

CALL  F  (XT  +  0. E¥HT ,Y ,P2 ) 

Y(I)=YT(I)+M*(.  7  07106781*  ( P2(  I )  -P  K  I)  )  ♦  P2  C  I) ) 

CALL  F (XT+FT  ,Y,P3) 

YCI)  =  YT(IMP1*C°0C  I)  **5  857  0643  8*Pl  (I)+3«41A21356*P2CI)*P3(  I)  )  /G.O 

GO  TO  X,  (6,9,11) 

RP  AX  =  0 
00  12  1=1, N 

RrPM^«n;!i?"-'«,"U«nIo,.,.3MY,,,-YPCin)/EK.x, 

RVAX=f MAXI  (F  ,FF\X) 

Y  (I)  =  Y  (I)  ♦  (Y  (I )- YP  (I)  )  /15«  0 
IFCRFAX-EFAX  )12, 13,17 
X  C=X  0  +  F 

IF(XO-X) 15 ,14 ,13 
RETURN 

IFIRFAX-G.  C  3 ’FFA  X)  30, 22,2 

H=HT 

XT  =  X0 

OC  19  I = 1 , N 
YF(I)=YT(I) 

YT(I)=YC(I) 

GO  TO  8 
ENO 
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SUBROUTINE  CCNIRl.  {  Y,  YP  ,T) 

OIMFNSICN  VC  1 2 ) ,  YP  Cl  2  > 

CCMMCN/CMA  IN/EECHD,AV,AVC5  » AVO  T  ,  J4  , '» 0  ,P  I  ,  TCPI ,  T  R  IMG  FT ,  NE  AR  , 
1YHPUS1 ,THRLS< , 41 , A  2, A3 ,E1* E2,3  3 ,T1 ,7  2,  FEC,  SIGHO  ,  THEC,  ANGLE » 
2GA,GB,GCtOE, CA,3R, TPT1 ,TRT  2,AL  t  PE, ALPHA7,£T4KTf CNEE  ,CLEE»V  CPT  » 
6SFE»rrE,VR,CVCR,  GNCR,CHCF,  THRU S,TRIM,EETAT,W71,HT2,XWT  ,CCNT,ZHT  » 
7AL0T, PEGT, ALCPT,  3E  C^T  ,  TCF-T  ,XT,  Y  T  ,?T,  IX  ,  I Y  ,  TZ  ,  IX  Z ,  I?  ,  HR  ,M  AS  S  ,  G  ,  S  , 
fiBfC»CX,CXCE  ,C2,CZQF,  CT,CK)E,CMQf  RVR2  ,ROH,STHE,CTHE,  TVGAIN 
REAL  IX, IY  ,I2,IXZ,  iR »  H  ASS 
CCUT =0. 

IF<T  «LT ,3. 1RETURN 

***  SPECIFY  CONTROL  HCCEt 

VOFT  =  C.  CONTROL  SURFACES  ONLY 
VO  FT  -  1 «  THRUST  VECTORING 
VO  FT  s  2 •  HOME  NT  THRUSTERS 

ANO  SFECICY  THRUST  CONTROL  GAIN.  TVGAIN=  %  CF  RLDOER  ANGLE 


**  SPECIFY  BETA  COKMANC  IN  DEGREES. 
TVGAIN=0.2 
VCFT=0. 

BECHO-C. 

,  CCNT=1. 

(  0E=-2E . 

DRHX=30. 

0AHX=15. 

8EC-CECMC/A V 
G1-. OS 
G2= .  5 
G3=. OB 

T 1-T  2= 1 C  C  C  C  « 

OE=DE-T+3. 

IFCOE.LT. -25.)  0S=-25. 

BEC=BECMC/ A  V 

DF=(G1*Y(6>-C2*C3E-PEC  >>  *A V 

DA=  CC3*Y  (<♦ ) )  #AV 

IF  CAPS  CCA) .GT.CAMX)nA  =  CA,0  A  MX/  APS ( OA ) 
IF (ABS 'OR) .G1  .CRMX )DR  =  CR*DRMX/ ARSC  CR> 
RETURN 
END 
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SUBROUTINE  V  E C TO  9  (  Y,  YF) 

OIHENFION  Y(1Z),Y«M12> 

CCMNCn/CNAIN/PFCMD,AV»AVC5,AVOT,JA  ,A  0,PI , TCPI ,  T  R IMOFT  ,  NE  AR  , 
lTHRUSl,THRl?o  31  ,A2,  A  2,B1,  e?,3  3  ,T1  ,T ? , F EC , SIGHC , THE C , A NGLE  , 
?GA;GP,GC»CE,CA  TRTi  ,TFT2t  AL  f  eE,  ALPHAT  ,STAKT,ChEE,CLEE»VCFT  , 
ESFEfCFE.VRfCYCRjCNCRfCFCPf THRU S ,TRIM , HE T AT , WT1 , NT2, XWT ,C CN  T , 2 F T  , 
7AL0T,BECT,aLCr-T,  BE  CP  T  ,  1  CFT  ,  XT,  YT.Z  T,  IX,  I V ,  IZ ,  IX  7 ,  IR  ,  WR  ,H  AS  S  ,  G  ,  S  , 
8B»C*CX»CXCE»C2»C?UE,C?SCNDE»CMQfRVR2fROH»STHE»CTHE»TVGAIN 
REAL  IX*IY,I2,IXZ, IR  j w  AS? 

IFIVOFT.EQ.O  .)  RETURN 
AlOEG=AL*A V 

IF(AIDEG,LT.  20.)  TVAU=C. 

IF(ALCEG.GT.  2C.)  TVAUs  < AlCE G-30 . ) /2 0. 

IF(ALPEG.GT.  EC.)  TVAU-1. 

IFC',  CFT.EC.2.)  GO  TC  99 
IF ( VCFT .£0 . 3  )  T VAU= 1 • 

B£TAT  =  T  VGA  IN  *T VAU*  CR/ A  V 
IF ( VOFT • EC . 1  . )  RETURN 
OR=  0  • 

RETURN 
99  CONTINUE 
ORMX=30. 

WT1=WT2=0. 
ang-tvgain*ccnx/av 
WTMX=(T1+TZ)  «\'T*SIN{  ANG)  /XWT 
HT=TVAU*WFX  »AcS  ( Q  F/  C  F  N  X) 

IF(DR. LT.fi.) FT  1=  WT 
IF(CR.GT.Q.) VT2=WT 
RETURN 
END 
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Appendix  E 

Aircraft  Modifications  For  Improved  Directional  Stability 

The  majority  of  the  high-performance  aircraft  built  to  date  has 
experienced  directional  stability  problems  which  failed  to  surface  prior 
to  aircraft  reaching  fully  operational  status.  As  previously  discussed, 
the  corrective  action  is  normally  to  impose  regulatory  limitations  on 
the  crew  in  order  to  avoid  that  portion  of  the  flight  envelope  vulner¬ 
able  to  departure  unless  the  severity  of  the  problem  necessitates  minor 
and/or  major  configuration  modifications.  The  following  are  examples 
of  modifications  to  the  airframe,  flight  control  system,  or  both. 

Airframe  Configuration  Changes 

One  of  the  first  problems  to  surface  stemmed  from  roll  inertial 
coupling.  The  corrective  modification  applied  to  the  F-86  and  F-100 
was  to  extend  the  vertical  tail  by  more  than  a  foot  in  length.  The 
same  technique  was  attempted  when  later  aircraft  developed  directional 
stability  problems  and  in  several  cases,  the  increased  tail  size  was  not 
only  ineffective,  but  proved  to  be  destabilizing  in  the  high  angle-of- 
attack  regime  (Ref  9  and  19) . 

There  has  been  considerable  research  conducted  to  determine  the 
effects  of  various  wing  modifications  on  directional  stability.  Drooped 
leading  edges  were  found  to  increase  both  directional  stability  and  lift 
characteristics  at  high  angle  of  attack  f  >r  certain  wing-body  combina¬ 
tions  (Ref  9) .  Although  drag  consideration  precluded  their  use  on  high 
subsonic  aircraft,  the  concept  led  to  the  retractable  leading  edge  slat 
such  as  that  used  in  the  F-4  Agile  Eagle  program.  Conversely,  there  are 
several  examples  in  Reference  19  of  configurations  using  leading  edge 
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slats  or  flaps  that  showed  little  improvement  in  directional  stability. 

Fixed  vertical  chin  canards  similar  to  those  used  for  Reynolds 
number  effect  correction  on  drop  models  have  been  used  on  a  few  aircraft. 
Their  effectiveness  was  limited  to  enhancing  spin  recovery  (Ref  3,  5,  9, 
21)  and  did  little  to  improve  directional  stability. 

Nose  strakes  such  as  those  to  be  used  on  the  Northrop  YF-17  light¬ 
weight  fighter  have  also  been  shown  to  improve  spin  recovery  capability 
on  model  tests  conducted  by  NASA  (Ref  10) .  This  improvement  was  only 
apparent  in  the  oscillatory  type  spin  and  no  effect  was  apparent  in  the 
flat  spin  mode.  Similarly,  no  evidence  of  improved  directional  stabil¬ 
ity  was  apparent  with  the  strakes  added. 

Control  System  Modif ieations 

Another  of  the  early  stability  problems  encountered  in  high- 
performance  aircraft  was  the  pitch-up  departure  of  the  F-101.  The 
pitch-up  occurred  prior  to  stall  and  was  followed  immediately  by 
departure  from  controlled  flight.  Obviously,  the  condition  was  a 
totally  unacceptable  risk  during  the  landing  configuration.  Ihe 
problem  was  initially  corrected  by  installing  a  stick-pusher  and  warn¬ 
ing  horn  (Ref  5).  The  stick-pusher  solved  the  problem;  but,  created 
new  problems  in  the  form  of  limited  maneuvering  capability  and  pilot 
dissatisfaction.  The  stick-pusher  was  later  replaced  by  the  Honeywell 
boundary  control  system  (Ref  5).  The  boundary  is  calculated  continuous¬ 
ly  by  the  system  as  a  function  of  Mach  number  and  compared  to  pitch 
rate,  angle-of-attack,  control  surface  position.  A  limit  function 
based  on  this  comparison  controls  the  engagement  of  the  system  which 
flics  the  aircraft  at  the  calculated  boundary  unless  the  pilot  over¬ 
rides  the  system  with  excessive  stick  force  or  reduces  the  stick  force 
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sufficiently  to  fly  off  the  boundary.  Unfortunately,  the  boundary  used 
is  that  defined  by  the  pitch-up  condition  which  was  never  corrected. 

The  result  is  an  artificial  maneuver  envelope  that  lies  well  within  the 
lift  limited  envelope.  The  region  between  the  two  envelopes  is  a 
safety  zone  that  without  the  stability  problem  could  be  vitally  needed' 
in  an  air  superiority  situation.  This  may  seem  an  unnecessary  require¬ 
ment  for  a  reconnaissance/air  defense  interceptor  aircraft;  but,  one 
must  realize,  the  weapon  system  is  worthless  as  an  interceptor  if  it 
cannot  defeat  the  aircraft  it  has  intercepted. 

A  current  example  of  a  control  system  fix  is  the  alpha-limiter/ 
beta-reducer  departure  inhibitor  system  to  be  installed  on  the  F-lll. 

The  alpha-limiter  senses  angle-of-attack  to  calculate  a  pitch  rate  for 
driving  the  pitch  damper,  to  reduce  command  augmentation  gain,  and  to 
increase  stick  force.  Tie  beta-reducer  feeds  differential  tail  as  a 
function  of  angle-of-attack  into  the  rudder  for  increased  roll  coord¬ 
ination  and  provides  yaw  rate  data  to  rudder  control  for  sideslip  reduc¬ 
tion.  The  name  of  the  system  is  the  best  statement  of  its  limitations. 
Like  all  control  system  cures  for  stability  problems  that  are  in  Ui.-e  to 
date,  the  alpha-limiter/beta-reducr.r  creates  a  buffer  or  margin  of 
safety  zone  at  the  performance  limits.  Obviously,  since  departures 
occur  in  critical  maneuvers  such  as  air-to-air  combat,  a  portion  of  the 
air  superiority  flight  regime  lies  in  the  buffer  zone.  The  loss  of 
this  portion  of  the  flight  envelope  may  seriously  hinder  mission  accom¬ 
plishment. 

Airframe/Control  System  Modifications 

As  mentioned  previously,  the  leading  edge  slat  configuration 
common  on  many  transport  aircraft  as  a  high  lift  device  has  been 
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successfully  used  &3  a  device  to  improve  directional  stability.  The 
Agile  Eagle  (F-4E)  modification  by  McDonnell  Aircraft  Company  resulted 
in  a  major  improvement  in  lateral-directional  stability  at  high  angles- 
of-attack  (Ref  5);  but „  ns  also  mentioned,  this  modification  has  not 
been  as  effective  on  other  configurations  (Ref  19). 

The  use  of  small  retractable  vertical  or  horizontal  canard  surfaces 
had  s4tiilar  effects  to  those  fixed  canards  mentioned  earlier.  While 
somewhat  effective  in  affecting  spin  recovery  (Ref  5  and  21),  they 
showed  little  evidence  of  improving  directional  stability.  Similar 
results  were  apparent  using  retractable  nose  strakes  (Ref  10). 

In  1959,  Cornell  Aeronautical  Laboratory,  under  contract  to  the 
Wavy  Bureau  of  Aeronautics,  successfully  improved  lateral-directional 
stability  of  the  F7U-3  by  installing  enormous  vertical  canard-mounted 
yaw  vanes  (Ref  1  and  5) . 

The  fact  that  strakes,  fences,  spoilers  slats,  canards,  oversized 
vertical  tails,  and  grotesque,  ungainly  yaw  vanes  are  required  to  provide 
some  semblance  of  directional  stability  in  various  modern  high-performance 
aircraft  is  sound  testimony  to  the  statement  made  by  Harold  Andrews  at 
the  1971  Stall/Post  Stall/Spin  Symposium  (Ref  5),  "We  still  do  not  have 
the  tools  to  design  the  a.'rcraft  right  in  the  first  place  ....". 
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